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NOMENCLATURE 




C^, thrust coefficient 

R rotor radius 

r radius to vortex 

induced velocity in wake in z direction 
X horizontal displacement 

z vertical displacement 

r circulation 

^ rotational speed 

c core size 

A6 defined in Fig. 1-2 
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SUMMARY 


This report is Volume III of a three volume series entitled 
**Free Wake Techniques for Rotor Aerodynamic Analysis** and 
coverin'j the following topics: 


Volume I (Ref. 3) ” Summary of Results and Background Theory 
reviews the results obtained to date using both complete 
and simplified wake models and siuimarizes the theoretical 
backgroxind on which these models are based. 


Volume II (Ref. 2) **Vortex Sheet Models** presents the 
results of computations using vortex sheets to model the 
wake and tests the sensitivity of the solutions to various 
assumptions used in the development of the models. The 
complete codings are included. 



Volume III (present volume) **Vortex Filament Models'* 
discusses results obtained using a vortex filament model, 
as opposed to sheets, again using various modelling 
techniques and including the computer codings. 


INTRODUCTION 


The application of computational fluid dynamics to the 
calculation of airloads on a hovering helicopter rotor is much 
more difficult than for an aircraft wing in steady flight. This 
is mainly because the trailing vortex geometry is more important 
to the in<3u:ed velocities on a helicopter blade than on a wing. 
Whereas a wing in straight flight leaves its wake far behind, a 
hovering rotor gathers it beneath itself in a spiral that remains 
close to the blade. Thus the induced velocities must be 
calculated everywhere in the wake as well as on the blade in 
order to give the correct bound circulation. The close proximity 
of concentrated vortices often makes convergence slow and the 
results extremely sensitive to vortex core size and other 
parameters. 

This report presents the results of airloads calculations for 
one case of a hpveri^ rotor for which good experimental data are 
available. The method used allows the trailing wake geometry to 
evolve freely according to the velocities calculated at discrete 
control points in the wake. This is known as a free-wake method, 
A loop calculates in turn the induced velocities, the geometry 
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and the bound circulation, iterating until convergence is 
reached. 

The method is loosely based on an original code written by 
J.D. Gohard (Hef. 1). That program has been developed by A. 
Tanuwidjaja (Ref. 2). In addition a simplified free-wake 
analysis has been done by R.H. Miller (Ref. 3 and 4) in which the 
wake immediately following the blade is made up of semi**inf initc 
vortex lines and the wake beneath the blade consists of circular 
vortex rings. In order to reproduce the experimentally observed 
bound circulation both these programs require assumptions in 
calculating the wake geometry and core size in the first 180 
degrees after the blade, as discussed in Ref. 4. This is a 
problem with the present method as well, as will be discussed in 
this report, and should become a focus of further research. 

The program presented here is a rewritten version of the 
Tanuwidjaja-Gohard program. In describing the wake, vortex 
sheets have been eliminated in favor of concentrated vortex 
filaments; some other minor changes have been made, and the 
program has been "cleaned up" to make it less confusing and to 
remove parts no longer used. 
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Section 1 of this report presents the free-wake method in 
detail. Section 2 discusses the results of the program for the 
tost case as different parameters are varied. The appendix 
contains the computer codes. 




OWG'Nft'- PRGE « 
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I. THE FREE-WAKE MODEL 

1* Method 

Each helicopter blade is treated as a concentrated line vortex 
of varying strength along the span. It is divided into discrete 
segments of constant vorticity where the strength of each segment 
is found from the angle of attack at its midpoint according to 
the formula: 

r • TTuac 

At the join of any two segments a trailing vortex is formed 
perpendicular to the blade with a strength equal to the 
difference in bound circulation on either side. (Figure 1-1). 
These trailing vortices are responsible for the downwash on the 
blade which, along with the blade pitch, give the angle of 
attack. 

Since the wake moves down as the rotor turns, the trailing 
wake forms a helix beneath the rotor plane. Each trailing spiral 
vortex in the helix is divided into straight vortex segments of 
equal angular length joined end-to-end. Thus the actual length 
of a segment depends on its distance from the hub (Figure 1-2), 
The control points where induced velocities are calculated are on 
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the vortex lines where two segments meet. The location of a 
particular segment in the wake is found by integrating the 
velocities up to that point on the vortex line to which it 
belongs. 


The whole wake is divided into three sections, the near, 
intermediate and far wakes (figure 1-3). The near wake retains 
all the individual trailing vortices from the blade to obtain the 
best resolution in calculating the bound circulation. After some 
azimuthal angle, usually less than 90 degrees, the intermediate 
wake rolls up the near wake into three concentrated vortices. 
The strength of a rolled-up vortex is the sum of the strengths of 
the near wake vortices contributing to it; its first azimuthal 
position is the centroid of the last azimuthal positions of the 
same near wake vortices: 


After several blade revolutions the intermediate wake is 
replaced by the far wake. The far wake is modeled as three 
semi-infinite vortex cylinders corresponding to the three 
vortices in the intermediate wake, (For a discussion of the 
effects of far wake modeling see Ref, 5, p. 16.) The start of 
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the far wake is one vertical spacing underneath the last 
azimuthal station of the intermediate wake, and the radius is 
assumed to stay constant* There are no control points in the far 
wake and no induced velocities are calculated there* 




The induced velocities due to the near and intermediate wakes 
are found using the Biot-Savart law for straight vortex segments 
with a correction for a viscous core size. The core allows the 
induced velocity to go to zero as a vortex line is approached 
instead of going to infinity* The induced velocity due to each 
segment in the wake and on the blade is calculated at every 
control point* These calculations consume the most computer 
time* A typical wake ^ contains two hundred segments or 40,000 
applications of the Biot-Savart relation for every iteration* 


The induced velocities due to the far wake are found from an 
elliptic integral series solution used by Miller (Ref* 3)* 

The self-induced velocity of a straight vortex segment is zero 
if a finite core size is used* But the trailing vortices are 
really locally curved, and a correction to take this curvature 
into account is needed* This has been derived by Scully (Ref* 5) 
based on the self-induced velocity of a vortex ring of finite 
core size found by Lamb (Ref, 6). The formula is: 
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= ^ (8-«n(|-tan(^)) - .25) 


where A6 is the angular length of the segments on either side of 
the control point. It is important to note that this is accurate 
only in cases where the arc length of the segments is large 
compared to their core size. This is true in the outer part of 
the wake where typical arc lengths are O.IH and core sizes are 
.OIR, but- it becomes questionable closer to the root. The exact 
geometry close to the rotor h»’b is much less important to the 
bound circulation, however, and these corrections are small 
compared to the total induced velocities. 

When overall convergence is finally reached, the main program 
calculates the lift coefficient. This coefficient, the final 
bound circulation and the wake geometry are compared with the 
experimental results. 


2. Parameters 

The important parameters to be varied in the program are the 
mesh sizes on the blade and in the near and intermediate wakes, 
the number of azimuthal stations in each wake, the criterion for 


roll-up in the intermediate wake, the tip loss factor and the 
viscous core size. 


a) Mesh Size 

The usual practice is for the bound vortex segments • n the 
blade to be long inboard and short near the tip where the 
circulation changes rapidly. In all the cases presented here, 
the distance between adjacent control points on the blade is .02R 
from the ninety-percent span out to the tip, .05R from the 
eighty-percent span to the ninety-percent span and .1 or .15R 
over the rest of the blade. In all, thirteen control points 
(twelve segments) are used. 

The near wake has the same number of trailing vortices as the 
number of bound segments requires; in this case, thirteen. The 
angular length of each segment is ten degrees. Usually the wake 
is carried 70 degrees after the blade. The intermediate vfake is 
made up of twenty-degree segments and is usually carried about 
740 degrees, or four blade passages, past the end of the near 


b) Roll-Up Criterion 


The tip vortex in the intermediate wake is rolled up from the 
outermost trailing vortex in the near wake to the point of 
maximum circulation on the blade. The middle vortex is rolled up 
three stations in from where the tip roll-up ends, or about the 
70% span. (As shown in Ref. 2 , the solution is not sensitive to 
the exact definition of this station.) The root vortex is rolled 
up over the rest of the blade, excluding the innermost vortex of 
the near wake. It is found that eliminating the root vortex in 
the near wake tends to smooth out the distribution of bound 
circulation on the inner half of the blade and makes convergence 
easier. 

c) Tip Loss • 

The tip loss factor compensates for the reduced lift at the 
very tip of the blade which is not well accounted for by the 
lifting line method. It simply sets the bound circulation at the 
last spanwise station equal to some fraction of the value it 
would have if the usual lifting line formula were used. This 
fraction is zero for all cases presented here, that is, the bound 
circulation of the last segment starting at 98% of the blade is 
set to zero; this is similar to the treatment of Ref. 7. 


d) Viscous Core Size 


The core size is probably the most important single 
It has a very large effect on the velocity Induced by vortex 
segments on nearby points. For example, if the core size of the 
tip vortex is * set equal to .02R, this causes the bound 
circulation to go up at the last spanwise station because the tip 
vortex no longer induces enough downwash there to pull it down. 
It makes no difference to the circulation over the rest of the 
blade. In this study, a core size of .OIR is always used at the 
tip; it may be larger in the rest of the wiahe. 

It is difficult to determine what core sizes' should be used in 
the program, and few experimental guides are available. In 
general, the results seem overly dependent on the core. One area 
that should be studied is how the "proper* core size depends on 
the mesh size. In the example given above, if the distance 
between stations were made larger, a larger core size could be 
used (and vice-versa) and yet the true core size should not 
change. For a further discussion of the effect of core size on 
the tip vortex in the near wake see Ref. 4. 



II. RESULTS 


Table 1 shows the characteristics of the rotor being studied. 
The experimental data were published in Ref. 8. 

a) Fundamental Solution 

The *£undamcntal solution" uses a core size of .OIR throughout 
the wake. It carries the near wake over 70 degrees and the 
intermediate wake over 740 degrees. The solution is given in 
Figure 2-1 along with' the experimental results; the bound 
circulation is plotted against the span and the tip vortex 
position is shown for every 180 degrees. 

There are three major problems with the solution. One is that 
the bound circulation is too low from about the 70% span out to 
the tip. This gives a lift coefficient (CT) of about -.00448 
instead of the experimental value of -.0046. Another problem is 
that the tip vortex remains too close to the blade in the first 
180 degrees. The third problem is that the wake defined by the 
tip vortex does not contract enough; it even begins to expand 
after 540 degrees. 




12 


b) Wake Contraction 


^ The wake does not contract enough either because the 
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intermediate wake produces too much radial velocity or because 
[the far wake produces too little. By removing the inner vortices 
in the intermediate and far wakes the wake contraction can be 
improved (Figure 2~2) but not without sacrificing bound 
circulation near the root, which yields a much lower lift 
coefficient. It is interesting to note that the circulation near 
the tip is unaffected. This implies that wake contraction has 
little to do with the rest of the problem. 


c) Core Size (Near Wake) 


The only way to raise the bound circulation near the tip is to 
increase the viscous core size in the near wake. All three parts 
of the wake, near, intermediate and far, contribute an important 
fraction of the downwash on the outer 25% of the blade. The far 
wake is too distant for any reasonable increase in the core size 
to have an effect. The core size of the intermediate wake can be 
important if it is larger than about .04R; this possibility, 
known as core burst, is discussed below. The near viake vortex 
segments are close enough for a small change in the core size to 
have a large effect. 


The result o£ increasing the core sise in tho near vaho (but 
notr as mentioned in Section X, at the tip) from ,01R to .02R is 
shown in Figure 2**3. The peak is higher and the tip vortox has 
accordingly moved down somewhat. The circulation in tho "valley” 
between the 70% and 90% span is still too low. 

d) Core Burst 

It has been observed that a concentrated tip vortex passing 
close to a blade undergoes a sudden expansion of its core (Ref. 
9). This fact can be exploited to improve the distribution of 
bound circulation. The intermediate wake tip vortex produces 
strong upwash on the blade fifom the 90% span out to the tip, and 
downwash everywhere else. Thus an increase in core size will at 
the same time raise the circulation in the valley and lower tho 
circulation at the peak. This is illustrated in Figure 2-4, 
where the assumed core burst is to .05R and the core size in the 
near wake is .02R. The wake geometry stays the same. 

In order to raise the peak back up again, the core size in the 
near wake must be increased to .05R. (It quickly reaches tho 
point whore a larger near-wake core size makes no difference.) 
This is shown in Figure 2-5, This last solution is the best 



obtainod 30 far. Tho lift coefficlont is *-.00453, about two 
porcont louer than tho exporimontal valuo. 
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o) Tip Vortex Location 

The doficioncy in circulation between the 70% and 90% span is 
possibly caused by the tip vortex coming too close to the blade 
at first encounter. The "natural” location of the tip vortex 
predicted by the program is closer to the blade than was observed 
in Ref. 8. It may be that the slow downward motion of the tip 
vortex is caused by inadequate modeling of tho roll**up in the 
near wake, possibly by not having enough trailing vortices at the 
tip. 

One way to tost this is to reduce tho near wake azimuth to 
only ten degrees - essentially rolling up tho tip vortex 
immediately. The results are shown in Figure 2-6 for the same 
case as Figure 2-5 (core burst to .05R and near wake core size 
.05R). Tho circulation and the lift coefficient are slightly 
larger, and the tip vortex location is a little farther from the 
blade. The solutions for other combinations of parameters are 
similar. This suggests that poor modeling of the near wake 
roll-up could be part of the problem. 
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In the earlier results from this program using a near wake 
spanning seventy degrees, the problem of the tip vortex location 
was much worse. It was discovered that this was because the 
program had not been allowed to converge properly. The initial 
approximation to the wake geometry does not roll up the vortices 
in the near wake. During the first few iterations, the main 
trailing vortex actually moves up. If the convergence test for 
the geometry is not tight enough the .program will converge before 
the near wake has had a chance to move back down. It usually 
takes over 60 iterations to converge within two percent. 





CONCLUSIONS 


The results of the program so far indicate that further 
investigation of the method is necessary. With a small core 
size I the bound circulation is too low on the outer part of 
fJie blade while the. tip vortex passes too close at first en* 
counter. The two problems are closely related. The slow 
descent of the tip vortex may be because of near wake modelling 
deficiv?ncy. Assuming a larger core size, corresponding to 
a burst vortex# results in better agreement with test data. 

A further area needing research is the relationship between the 
core size and the mesh size. 
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Table 1. ROTOR CHARACTERISTICS 


Straight 

Blades - 2 
Twist - 11® 

Collective pitch at 75% Radius 9.8® 
Solidity .0464 






















'jiF^ENDIX 

LIST OF VARIABLES 


Follouing is a list of important variables in the FreeWake 
program. 


alp 

alphas 

blades 

cd 

coef £ 
dpsii 
dpsin 
etanv 


angle of attack, alpaflamda-thetac 
angle of attack at stall 
real number of blades 

(cdO,cdk) drag coefficients, cd«»cdO + cdk*alp*alp 

(coeffl) tip loss coefficient 

angular length of intermediate wake segments 

angular length of near wake segments 

spanwise location of endpoints of bound segments 


etaiv 

etal 

cta2 

eta3 

facgam 

facgeom 

facvel 

f lamda 

fmu 


initial spanwise location of intermediate 
wake trailing vortices 

spanwise location of intermediate wake trailers 

relaxation factor for LOOPl 

relaxation factor for INTGR 

relaxation factor for L00P2 

inflow angle due to downwash on blade 

net inflow — used only for wind turbine cases 
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o 


fpsl ' 

£ps2 

£ps3 

game 

gamt 

gamti 

gamt ip 


core size £or near wake 

core size £or tip vortex in near wake 

core size o£ intermediate wake (core burst) 

bound circulation 

trailing vortex strengths in near wake 
rolled~up vortex strengths in intermediate wake 
same as gamti (3) 


int 

itest 

iwrite 


knnvr 

knivr 

liml 

lim2 

IP 

It 

nblds 

ncase 

niter 

niva 

nivr 


parameter for intermediate output: int»l prints 
li£t coefficient and wake geometry for every fifth 
iteration 

itest«0 if all subroutines have converged 

parameter for output: iwrite»0 gives only final 
results, iwrite»l gives results for every iteration, 
iwrite»2 gives induced velocities from each part of 
wake (near, intermediate and far) for each iteration 

initial value of nnvr 

initial value of nivr 

maximum number of iterations for LOOPl 

maximum number of iterations for main loop 

power coefficient, or drag integrated over span 

lift coefficient, or lift integrated over span 

(nbldsl) integer number of blades 

number of different cases in file to be run 

current iteration number of main loop 

number of azimuthal stations in intermediate wake 

number of spanwise stations in intermediate wake 






c 
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nnva number of azimuthal stations in near wake 

nnvr number of spanwiso stations in near wako 

sigma solidity of rotor 

theta blade pitch at endpoints of bound segments 

thetac blade pitch at midpoints of bound segments 


wxiv x,y,z components of induced velocities in 

wyiv intermediate wake 

wziv 


wxnv x,y,z components of induced velocities in near 

wynv wake 

wznv 

xiv x,y,z coordinates of control points in 

yiv intermediate wake 

ziv 

xnv x,y,z coordinates of control points in near wake 

ynv 

znv 


At the end of the computer code a sample output is given (for 
Figure 2-5) showing the lift coefficient, bound circulation, 
induced velocity on the blade and wake geometry. 

The bound circulation (GAHC) and induced velocities (WZC, W7C, 
WZC) arc given at the centers of the bound vortex segments 
whereas the trailing vortices (which give the wake geometry) are 
located at the ends. 

The integer headings in the wake geometry refer to the 
azimuthal station number. The first line under each heading is 
the azimuthal angle (PSI), the second is the radius (R) and the 
third is the vertical distance (Z) from the rotor plane. The 
first azimuthal station is at the blade, so here PSI»0., Z^O. and 
R gives the spanwise positions of the ends of the bound vortex 
segments. 
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ORIGINAL PAGE B 

OF POOR QUALiri 


Computer Codes 


FREE-WAKE ANALYSIS OF A WIND TURBINE AERODYNAMICS 


common/pana/ Iwr . Ird, Hraa, Hml , 1 in<2;ni ter. Iwri te 
coremon/Qooni/nblds 1 * nbl ds« s Iona, ff«u, oton( 15 ) * Wnnvr . ota 1 ( 6 1 * kni vr • 

S 1 tw< St , thotadi 15 ) , theta( 15 ) . thotac( 14 ) , thatoO. thotOd. a Iphns . 
ScdO.cdk.dpaInd.ctpsIn.dpsi Id.dpsi 1 .coof f .coof f 1 ,c.s. blades, 

Snnvr , nnva , nnvr 1 , nnva 1 , nncr , nnca,atanv( IS) , otanc( 16 ) , 

Snlvr.nlva.nlvpl .nival ,n1cr,nlca.atalv(6).etalc(7). 

Sntva,ntva1,ntca.fps1,fp32 „ 

coniaion/samma/oamc( 16 ) . gamtl 15 ) . gamt 1 ( 3) .gamt Ip . eta I .ata2. ota3 . 

&k1.k2,k3.facga.'9 , ..... 

cotitflion/posl t/xnv( 15. 18) ,ynv( 15. 13) ,znv( IS. 18),xlv(6.S0).ylv(6.50). 
&z1v(6.S0) 

coti«iion/va1oc/wxnv( 15. 18) ,wynv( 13. 18) .wznv( 15,18).wxlv(6.50). 

&wylv{6.50).wzlv(6.30),facgaon 
coimron/ Ini t 1al/wxc( 14) , wyc( 14 ) ,wzc( 14) 

cotKi»on/v I rvdo/ If ar. nf ar.xl .x2,y1.y2,z1.z2.x.y .z. facval .gm.ux.uy .uz 
comnon/aot f / 1 sol f 

common/ 1oopsave/wxv( 16) .wyv( 16) ,wzv( 16) 
common/spoc I a 1 / I tgr 
com«on/extraspoc/tps3 

name) Ist/fparm/ Iwr . Ird. 1 Ins, 1 1ral. 1 Ib 2 . Iwrita.nfar. Ifcr.gam.wx.wy. 

Swz.facgam. facval , facgeora.ncaso. Isel f . Int, I tgr 
real I Ip. Ip.- It 
data p1/3. 14159/ . 

Iwr«6 

1rd»5 

lself-0 

coeff*0. 

nb1ds1«2 

alphas*!. 

cooff 1".S 

11ffl8*S0 
I 1ml*SO 
1lm2-15 
nfar*90 
1 f ar«0 
lnt“0 

read! Ird.fparn) 
do 901 nca 5 » 1 ,ncaso 
call inputl(ncas) 

lf(fps3.no.fps2) wr1to( lwr.53) fps3 
format!" CORE BURST TO ".f5.2) 

If ( Ifar.eq. 1 ) write! Iwr, 54) 
format!" ELLIPTIC INTEGRAL FAR WAKE ") 
write! Iwr. 57) cooffi 
format!" TIP COEFFICIENT ".fS.2) 

If! Itgr.aq. 1) write! Iwr, 58) 

format! “ ROLLED UP TRAILERS MEET AT POINT") 

do 10 l»1,nnvr 

et*etanv! l ) 

wxc! 1 )«wx*et*ot 

wyc!l)“0. 

wzc! I )"ot«wz/.8S 

game! 1-*-1 )«gam»Qt/.9l 

If !ot.gt. .91 ) wzc! 1 )*lwz*et/.83)»! 1 , -et)/. 15 
If ! Qt .gt . .91 ) game! I-*-! )■( gam»et/ . 9 1 ) •! 1 . -et)/.09 
cont Inua 

wzc! nnvr )»wzc! nnvr- 1 ) 



I 


^'v •: 


^RKxIkml i5 

OF POOa QUALiTV 
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c • 

c • 

c • 

c • 




t004 

1 

2 

c 


150 


160 

1003 


gamc(nnvr)*0. 
ganic( t )«0. 
g£rnc(nncr)«0. 
call compa 

MAIN LOOP 


Vwplta»lwrlto 
do 1 n1tsr*1,11ffi2 
I toat*0 
facgoona.S 
facvsl B.O 

If (niter. to. 2) facgsora«1. 
if (kwrIte.cq.O) Iwrito^O 
f 1 tor "nl tor/5, 
fnl tor»float(nl tor/5) 

If ( fnlter.oq.fi tor.and.kwrIto.oq.O) lwrlte«l 
If (nl tor.eq. 1 .and.kwrlto.eq.O) Ivirlte-I 
call loop2(ktest) 

1 toat>ktost 
call intgr(ktost) 

I tost* I tost«^kte9t 
call loopi(ktest) 

1 teat* 1 tost+k test 

If ( Itest.eq.O) goto 2 

1 f ( int .oq. 1 .and. Iwr 1 to.gt .0) goto 1003 

continue 

continue 

continue 


wrlto( Iwr, 100) 
wr 1 te( Iwr, lOB ) 
wr 1 te( Iwr, 120) 
wrlte( Iwr, 120) 
wrl te( Iwr, 120) 
wrl to( Iwr, 120) 
wrl to( Iwr, 107 ) 
do 150 J ■ 1 . nnva 
wr 1 te( Iwr . 109 ) 
wr 1 te( Iwr, 120) 
wrl te( iwr . 120) 
wr 1 te( Iwr, 120) 
cont Inua 
do 160 J"1.nlva 
wrlta( Iwr, 109) 
wrl to( Iwr. 120) 
wrl te( Iwr. 120) 
wrlto( Iwr, 120) 
cont Inuo 
cont I nua 
lt*0. 


1 teat 

(gamc( 1 ) , I*! ,nnvr) 
( wxv( I ) , 1 * 1 . nnvr ) 
(wyv( 1 ) , 1*1 .nnvr) 
(wzv( I ) , 1*1 .nnvr) 


) 

(xnv( 1 , J ) . I"1 .nnvr) 
(ynv( 1 . J ) , l"1 ,nnvr) 
(znv( 1 , J ) , I"! .nnvr) 


J 

(xlv( I, J ). l"1,nlvr) 
(yiv( 1 , j ) , 1-1 ,nlvr) 
(zlv( I . J ) , 1« 1 ,nlvr) 


lp«0. 

do 148 l«1.nnvr-l 

et".5*(otanv(l)+Qtanv( 1+1)) 

wi«wzv( 1 )+fnu 

wy"wyv( 1 )*at 

f 1 amda«atan( wz/wy ) 

u»sqrt(wz**2+wy»*2) 

a1p«f 1 amda- thetacl l ) 
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33 


I 


148 


901 

100 

109 

108 

106 

107 

120 

999 


CdBCd0«cdk«alp*a2 

If (AOS (alp) .gt .fids (alphas ) } cd>2. •cdO«^cdk«a1p«*2 

I lp«u*ganc( 1^ i ) 

tHp«Hp*cos(f lareda) 

f1 lpa*l lp«sln(f latnda ) 

dp>u«u*pl •( stgnia/&1adas)*cd« .5 

fdp*dp*CQS( f larada) 

top"dp*sln( f Ifiinda) 

tp*tl Ip+tdp 

fp»f 1 ipf'fdp 

I t*l t>tp»(atanv( 1*1 )«otanv( I ) ) 

1pa1p<f fp«(otanv( 1 + 1 )»"2-etanv( 1 )*»2) 

continue 

1 t*l t*bladas/pt 

lp«lp*dlfides/2./p1 

wr1ta( Iwr. 108) It, Ip 

1 f( Int.eq. 1 .and. I test. na.O. and. niter. 1 1. 1 11x2) goto 1004 
lwrite"2 

call 1aop2(ktest) 
continue 

format(//.' MAIN FINAL RESULTS -- CONVERGENCE 14 ) 
rormat(/.14) 

fortnat(/,* GAMC, WXV. WYV.W2V: * ) 

format(/.“ LT--.ft0.6." LP«".f10.6) 

forraat(//,* WAKE GEOMETRY:*) 

forraat(2(9f 10.6./)) 

cont Inue 

stop 

end 











rf!J/U/C> *K.' 


c 

11 

77 

c 


65 

66 
12 
13 


suopoutina Input Kneas) 

common/parm/iwr . Ird. I tma, 1 1m 1 , 1 1m2. nt ter , Iwp 1 1* 
COfflfflon/gQpm/nbldsI .nbtds.sfgma.fmu.«tan( 15) .knnvr .atal (6) . knJ vr , 
, & I twist , tbptadC 15). that a( IS) , tnetac( u ) , thetaO. thotOd. alphas. 

&cd0.cdk .dpslnd.dps1n.dpsi td.dpsi i .coaf f .coeff 1 .c.s. blades. 
6nnwp.nnva.nnvp 1 .nnval , nncp,nnca,etanv( 15 ) .etancl 16 ) , 

&nlvp. nWa. nivrl , nival . ntcr.n1ca.ota 1v{6) .eta1c(7) . 

Sntva.ntvs 1 .ntca. fps 1 . fp32 
corareon/extraspec/fps3 

equivotenco (nnvr.nncrl ) . (nlvr.nlcrl ) . (nnva.nncal ) 
equivalence (niva.nicat ) , (nnvrl .nncr2) . (nivpl .nlcr2) 
name I Is t /case/nb 1 da . s 1 gma . knnvr . nnva . dps 1 nd . e t an . kn 1 vr . 
&nlva.dP3i id.etal , fmu. thetad. alphas. CdO.cdk, 

&fps1 ,fps2.coaff 1 .fps3 
data conv. twopi/. 017453293.6.263185308/ 

I 

read( Ird. case. epp«999.aiid«aaa) 

dps ln«dpa lnd»conv 

dpsi iadpsi id*conv 

do 11 1"1, knnvr 

theta( 1 )athatad( 1 )*conv 

continue 

b1ades«f loat(nblds) 
c"Cos( twopi/blades) 
s*s1n( twopi /blades) 
continue 

nnvr»knnvr 
n1vp«knivp 
nncp*nnvp+ 1 
nlca^nlva+l 
n 1cp»nlvp*1 
niva 1 "nlva- 1 
nnvp 1 -nnwr- 1 
n1vp1»n1vp- 1 
do 65 I • t ,nnvr 
etanv( 1 )«etan( 1 ) 
cent Inue 
do 60 I■1.n1vp 
etalv( 1 )-Qtat( 1 ) 
continue 
do 12 l«2.nnvr 

etanc( 1 }■ (etanv( 1 l^etanvC 1 - I ) ) * .5 

continue 

do 13 l«2,n1vp 

otatc( 1 )*(eta1v(1)'»etalv( 1-1))*. 5 
cont Inue 

otanc( 1 )• 1 ,3*etanv( 1 )-.5*etanv{2) 
eta1c(t)"l.5*ota1v(l)-,5*etalv(2) 
etanc( nner )■ 1 . 3*etonv( nnvr ) - . 5*otanv( nnvr- 1 ) 
otalc(nlcp)*1 .5*etaiv(nlvr )-,3*etaW(nivp- 1 ) 

nnva1»nnva-l 
nnea^nnva* 1 
do 32 I*t.nncp2 

thetad 1 )*theta( i )*( theta( I ♦I )-theta( ♦ ) )* 

4(etand l••■1 ) - e tanvl 1 ) )/( etanv( l■••1)-etanv(1)) 
cont Inue 

1 f (noi ds . 1 t . 1 . op.nblds .gt . 8 ) goto 60 
If (nnvr. gt. IS. or. nnvr, It. 2) goto 60 



32 


35 




if(n<vr at.6.or.nwr.1t.2) goto 60 

<f(nnvB.gt. 18 .or.nnva.Jt. 2 ) goto 60 

t f (n« va.gt .50-or .n W« . 1 t . 2 ) g 
«f(ctanv( n.lt.O. ) goto 60 

do 80 <«a.nnvr 

lf(atonv( l).le.atanv( «-D) goto 60 

continue 

i?(JUw?;)Ii;.etnW(.-0) goto 60 
continue 

& *• MU *fH*5) * 

ft/ •jK.'lf 10.5) ^ 

wr1te( «wr. I53)a1pna».cd0.cdk 4 -4.- . fS. 3. ••ALPHA**2* . 

fnrmfltC STALL ANGLE ' . M 2 . 5 . / . * CO" = 

wr1te( twr. ’54 ) Urns. Hnl . j inj. fpsi . fpa 5 E„pi .LOOP 1 .L00P2: • .3 14, / , 

wr 1 te( iwr . 1 55 )knnvr . nnva . ( *tan( ‘ ) • . . , j. - ) . • . Af 10.5 . 

fornatC* NEAR WAKE DEFINITION. ( .12. 

ft/,5«.6M0.5./.4t 10.5) I.) knivr) 

wrltel1wp.156)kn»vr.ntva.(etal(1).l ^ -).-.4f10.5. 

^ornate INT. WAKE DEFIN TION: C . 12. . .<2. ). 

ft/. 5 x. 6 f 10.5./.5x.4f 10.5) 

return 

-“Input ... e«"Or on liRO) •) 

Stop 


TNOOf riLE ON (IR01-) 

Stop 

fornat( 8 f 10 . 5 ) 
f ormatl 312) 

continue 

rorJ.lt! INPUT DATA INVALID OR OUT OF RANGE *) 

wri to( iwr .case) 

stop 






PR<5£® «■•• S'vj 

QU».rrf ’'•■'I V. , : ' 

' ■ - ■'» ;.•) 

of 

\ 



subroutine conpa 

comnon/parm/iwr. ird. f 1ms. 1 Imt, I Ira2,nttap. Iwrlto 
common/geom/nbldsl .nblds.slona, fmu,atan( 15 ) ,knnyr,*ts< 1 C ) .-u-r 




4ltwlst,thetad( l5).thota( l5),tnotacf 14). thataO, thatOp.a un<jL 
&cdO,cdH.dpslnd.dpsln.dps1 id.dpsi 1 .coaff .coa^f 1,c.s.piaa«»t> . 




Snnvr.nnva.nnvri .nnval .nncr.nnca.otanv( 15) .otsnc( 16). 
&ni vr.ntv8.nl vri .nival .nicr.nica, eta 1v(6) .eta ic(7) . 


• 


Antva.ntva 1 ,ntca. fps 1 . fps2 

common/post tAnv( 15. I8).ynv( IS. 18) .znv( IS. 18) .xlvfe.&O) . 

&y lv(6.S0) .z tv(6.50) 
common/ Ini t 1a1/wxc( 14 ) . wyc( 14) .wzc( 14) 
dimension wrn( 15 ) .wr i (6 ) . wzn( iS) . wz 1 ' 5 ) 
equivalence (nnvr ,nncr 1 ) . (n1 vr.nicr 1 ) . (nnva.nncal > 
equivalence (nlva.nica) ) . (nnvrl .nncr2) . (nlvrl.nlcr2) ' 
data pl/3. 14159/ 



e 

• 




do 21 l« 1 .nnvr 
xnv( 1 , 1 )«etanv( 1 > 
ynv( 1 . 1 )«0. 
znv( 1 . 1 )«0. 

pni«o. 

r»etanv( 1 ) 

do 21 J"1.nnva*1 

pni l»phl+dp3 In 

r 1 "r+dps 1n»wxc( 1 ) 

xnv( 1 . J+ 1 )"r l*cos(pnl 1 ) 

ynv( 1.J+t)»ri*s1n(pnli) 




znv( 1 . J + t )»znv( 1 . J )«'dps 1 n«wzc( 1 ) 
phi •phi 1 




r»r 1 



21 

cont Inue 




1 f ( iwr 1 to. 1 1 .2) goto 103 
wr 1 to( Iwr. 500) 



500 

format(//. INITIAL WAKE GEOMETRY:*) 
do 9 1 J • 1 . nnva 
wr 1 te( Iwr .505 ) J 

»r 1 to( Iwr . 5 10) (xnv( 1 . J ) . 1 ■ 1 .nnvr ) 
wr 1 te( Iwr ,5i0) (ynv( 1 . J ) , t«1 .nnvr) 
wrtto( Iwr. 5 10) (znw( 1 . J ). l«l.nnvr) 



91 

continue 



103 

continue 



505 

format(/. 14) 



510 

format Mx,2(5x.9f 10.6./) ) 



c 

• 



c 

• NODES AND VELOCITIES FOR THE INTERMEDIATE WAKE 



c 

• 




do 30 1*1 . nivr 

et*etaiv(l) 

do 31 I2«2.nnvr 

I f (ot . 1 t .etanv( 12) ) goto 32 



31 

cont inue 



32 

continue 

a»(etanv( 12) -at )/(ctanv( l2)*otanv( 12*1)1 
wz 1 ( 1 )»a*wzc( 12- 1 )♦( 1 . -a) »wzc( 12 ) 
wr 1 ( 1 )«a»wxc( 12-1 )♦( 1.-a) »wxc( 1 2 ) 
x1v( 1.1)«a*xnv( 12-1 ,nnva)^( l . -a)*xnv( 12. nnva) 
ylv( 1. 1)»a»ynv( 12-1. nnva )+( 1,-a)»ynv(12. nnva / 
z 1 v( 1 . 1 ) "a-znvl 12 - 1 .nnva )♦( 1 . -a ) *znv( 1 2 .nnva) 
ps 1 "dps in»( nnva- i ) 
r«3qrt(x1v( 1 . 1)«*2*ylv(1, 1)*»2) 









oo 22 J»1,rt!v>a«l 

psH-pst+dpsn 

wrr«wri(<) 

If (ps< 1 .gt.pl ) wpr«wrp/3, 

I f (ps 1 1 .gt .2. *pi ) wpp«wpp/2. 

I f I ps I 1 .g t . 3 . *pl ) wPP«wpp/2. . 

I f (ps I I .gt .4 . *p« } wpp*wpp/2. 

Pl»P*CP8l l»WPP 

* »v( I . J ♦ t )»P I *C09(pB 11 ) 
ylvl I . 1 )ar 1 csinlps 1 1 ) 

wz«wz1 ( I ) 

If(ps1 l.gt.pl . and, I . eq , n 1 vp ) wz*2 . *wz 
ziv( 1 , j«i ) bz1v( 1 , j )>wz*apsi 1 
pap 1 

pa laps 1 1 
22 continua 

•30 continua 

lf( Iwpito. It.2) goto 104 
oo 92 J«l,nlva 
WP1 ta( lwp.505) J 

wpitadwp.SlO) (xlvd.jj.fi.nlvp) 
wpl tadwp.SlO) (ylvd , j ) . |at .nivp) 

^ — wpltadwp.aiO) (zlvd.j).iai,nivp) 

■ "2 conttnuo 

104 continua 

patupn 
and 


o 
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subrouttn* )oop2(kteBt) 

Thl* subrout ina calculates the Inauceb valocittes at all 
tba points in tha near and Intarriaalats wake. Only tralllno 
signent elements are used. The wake Is rollod up Into three 
traiM^ vortices starting at tha intermediate wake. Only the 
ends OP the segments are considarod -- the positions of the 
centers are never used or even calculated. 

real k3 

COmmon/parm/lwr, Ird, I Iras, I let, I In2, niter, Iwrita 
common/geom/nbldsl.nblds.slgma, fmu.otan( 19 ) , knnvr .etal (6) .knivr. 
el twist , the tad ( 15) , theta( 15) , thatacl M ), thetaO. thetOd, alphas. 
acdOtCdk.ops ind.dps in.dpst id.cpsi i .coef f. coePf I ,c, s.blades. 
annvr . nnva . nnvr I , nnva 1 , nner . nnca , etanv( 15). e t anc( 1 6 ) , n I vr , 
aniva.nivr I .nival ,nlcr,nlca.otalv(6),etaic(7), 

Antva.ntval .ntca. fpsi ,f?s 2 

common/ gamma /game ( 16) .garnt ( 19 ) .gant I (3) .gant Ip.etaO) . 

6 j 1 .kkl 2 ) . faegam 

common/ve loc/wxnv( 19. I8).wynv( 19. 18) .w*nv( 19. 18 ) .ws lv(6.90) . wy lv(6.50) . 
aws tvt 6 , SO) , Paegeom 

common/vlnoo/IPar.nfar.x1,x2.y1.y2.2l.22.x.y.x.Pacvol.gm.UK.uy.u* 

comnon/post t/xnv( 15,10) ,ynv( 15, t8),znv( 15 , 18 ) « x1 v(6 ,50) ,y I v( 6, 50) , 

52 1 V( 5 , 50 ) 

common/save/wxnvs C 1 9 . 18 ) . wynvs (19.18), wznvs ( 1 3 . 1 8 ) . wx | va ( 6 . SO ) . 

&wy tvs(6.50) .w*lvs(6.90) 
common/ael P/lselP 

CQmman/oxtrasoec/Pps3 
data twopl . fpl/6.2Q3ia. .079377/ 

TRAXUNG VORTEX STRENGTHS CALCULATED BY SUBROUTINE ROLLUP 

call rollup(2) 

IP( iwrite. It. Dgoto 3 

write! Iwr, 100) niter 

wrlte( Iwr, 123) (gamc( I ) . I»1 .nnvr) 

wrlte( Iwr. 123) (gamt( i ) . l-l .nnvr) 

wri te( Iwr, 123) (garntK I ). l-l.nivr) 

wri te( Iwr, 1 10) (eta( I ) , I«1 .3) 

continue 

initialize arrays for new inouceo velocities 

do 10 1»l,nnvr 
do 10 J»1,nnva 

wxnv(l.))»o. 
wynv( I , j )»o. 
wznv( i , j ) "O, 

cont I nue • 

do 12 l•t,nlvr 

do 12 J«i,niva 

wxlv( I . J )«0. 

vyiv(l,j)«0. 

wr iv{ 1 , j )«0. 

continue 

velocities inouceo by blaoe bound circulation 

X 1 ■etanv( I ) 
y 1 "0. 

21"0. 





ya»o. 

ia«o. 

do 30 1«t,nnvr-1 
»3*«t«nv( ) 

Qffl*«ganc( ) 

ni»* 

n2«t4i 

call vindn(nt .n2.2) 
call v<ndt(nl,n3) 

X i«m2 

20 ccntinua 

(f( 1wrlto.no. 2) goto 32 
do 21 ]«l,nnva 
wri to( iwr , 1 19) J 

wr1to( Iwr, 120) lwxnv(k, J ) ,K» 1 ,nnvr ) 
wr 1 to( Iwr , 120) ( wynv(k. J ) .k» 1 .nnvr ) 
wr I to( iwr, 120) (wtnv(k, j ) ,k» 1 .nnvr ) 

21 continuo 

do 23 J"1,n1va 
wrlto( Iwr, 1 19) J 

wri te( Iwr, 120) ( wx1v( 1 , J ) , 1 «1 , nl vr ) 
wr 1 te( Iwr , 120) (wy1v( 1 , J ) , 1"1 ,n1vr) 
wr 1 to( Iwr , 120) (wz1v( 1 , J ) , 1 • 1 , nivr ) 
23 continuo 

22 continue 
C 

C VECQCITIES INOUCEO BY NEAR WAKE 

C 

dps2*fps2 
do 30 I ■ 1 ,nnvr 
fps2*fpsi 

1 f ( 1 .oq.nnvr- 1 ) rps2«dps2 
do 30 J«1,nnwo-1 
Kl»xnv( 1 , J ) 
y1«ynv( 1 , J ) 

X1«xnv( 1 , J ) 
x2«Mnv( 1 , J*1 ) 
y2«ynv( 1 , 1 ) 

X2«znv( 1 , J+1 ) 

gmagantl I ) 

nl ■{ J • 1 )«nnvr+1 

n2*nl+nnwr 

call vlndn(nt ,n2,2) 

1 f I 1 .oq.nnvB ) n2»0 
call vtndl(n1,n2) 

30 continue 
fps2*dos2 

1 f ( Iwr I to.no. 2 ) goto 32 
do 31 J» 1 ,nnvB 
wri te( iwr, 1 19) J 

wr 1 to( Iwr , 120) ( wxnv(K. J ) ,k»1 ,nnvr ) 
wr1to(1wr,120) ( wynv(k . J ) ,k« 1 ,nnvr ) 
writal Iwr, 120) (wznvlk.j ),k*1, nnvr ) 

31 continuo 

do 33 j*1.nlva 
wri ta( Iwr, 1 19 ) J 

wr1t®( Iwr. 120) (wx1v(1,J), 1«1,nlvr) 
writQllwr,l20) IwyW(I.J), 1-1, nivr) 
write! Iwr, 120) (wzlv(l.j), I•1.n1vr) 
33 continue 



40 






continue 

INTERMEDIATE WAKE 

ph1b"twopt /blades 
do 40 l■ 1 .n^vr 

tf(t.aq.l) ksaveaniva 
lf(t.eq.2) Ksave>nwa 
1f(1.eq.3) kSBva«n(va 
do 40 J*i.ksavo*t 
xt«xtv( 1 . J ) 

Vi«y lv( I , J ) 
zl«xiv( I . J ) 
x3*x w( 1 , ]♦ t ) 
y2»y 1v1 1 . O 
Z2«x1v( 1 . j^l > 
gm"gamt1(l) 
dp82*fps2 

ph1l"f loat(nnva-4 )»dpsin ♦ f loat( J -1 )*dpa1 1 

If (pni t .gt.pniD-dpsIn) fps2«fps3 

nl*( J- 1 >*n1vr nnva*nnvr + 1 

n2»n l^niwr 

If(J.eq.i) n1«0 

call v1ndn(nl .n2.2) 

call vlndl(nl.n2) 

fps3*dps2 

continue 

If ( lwrite.ne.2) goto 42 
do 41 J • 1 , nnva 
wri to( iwr, 1 19) J 

UP 1 to( iwr , 120) (wxnv( 1 , J ) . !■ 1 ,nnvr) 
wr1t0(1wr,12O) ( uynv ( 1 , J ) . 1 ■ l . nnvp ) 
wp I te( Iwr. 120) (wznv( 1 . J ) . 1 • 1 ,nnvr) 
continue 
do 43 j»l,nlva 
WP 1 ta( iwr, 1 19 ) j 

wr 1 ta( twp, 120) { wxiv( I . J ) , 1 • 1 ,n1vr) 
wr1ts( Iwr, 120) (wylw( 1,j ), 1"1,n1vr) 
wPl te( iwr, 120) (wzW( i , j ) , 1"1 .nlwr) 
continue 
continue 

FAR WAKE 


1 f ( 1 f ar.eq. 1 ) goto 45 

dz«zlv(nlvr,nlva)-zlv(n1vr,n1va-1> 

gm«gamt1p 

x1"xlv(nlvr,n1va) 

y1“ylv(nlvr,nlva) 

Zl*z1v(nlvr,n1va) 
phi 1 >atan2( y 1 , x I ) 
r«sqrt(x1*x1 ♦ yl»y1 ) 
do. 50 J • 1 , nf ar 
phi l»phl 1+opsl 1 
X J«r»CQs(pn1 I ) 
y2"r*s ln( phi 1 ) 
z2*z l^dz 

nl«J ♦ nnwa»nnvr ♦ niva«n1vr 
n2«nl+1 

call vlndn(n1 ,n2,2) 


41 





catt v(nd«(nt,n2) 

Kt«x2 

Vl»y2 

zi»z2 

90 continuQ 
goto 49 

<9 continua 

fkapa*phlb/dps( i 
kapa«(Mx(fkapa * ,S) 
do 300 k1»».nlvp 
Kl*xlv(kl .niva) 
y lay lv(k( ,niva) 

H-zlv(kJ.nWa)+( 2 iv(ki,ntva}- 2 lv(k«.niva-kapa)) 
aganagamt t(kl)/(2tv(tci.ntva)-zlv(k( ,n(va-kapa) ) 
ra(5qrt(xt«xl ♦ y1*y1)) 
do 46 I a 1 , nnvr 
do 46 jai.nnva 
na t 

xaxnv(l.j) 
yaynv(».j) 
zaznvfl,}) 
z3aabs(2l - 2 ) 
phlaatan2(y.x) 
ata3asqpt(x«x ♦ yay) 
goto 200 
210 continuQ 

wxnv( 1 , j )-wxnv( 1 . J ) + wr«cos(pb 1 ) 
wynv{ 1 ,j)pwynv(i,j ) ♦ wraafn(phl) 
wznv( f , J )pwznv( 1 , J } ♦ W 2 
46 continue 

do 220 lai.nlvr 

If(l.eq.t) kaavoaniva 
lf(1.eq.2) ksavaaniva 
1^(1. eq. 3) ksavaaniva 
do 220 Jal.ksava 
na2 

“•xlvd.J) 
yaylvf I ,j ) 
z-zlvd.j) 
z3aabs(2i-z) 
pni-atan2(y,x) 
ata3*sqrt( x«x ♦ y*y) 
goto 200 
230 continua 

wxlv( 1 . j )awx1v( 1 . J ) wr*cos(pht) 
*^y«vd.J)-wy1v(1.J) ♦ wr-sln(phl) 

WZiv( 1 , J )av,zlv( I , J ) W2 

220 continue 

If ( twrtta.ne.2) goto 32 
do 51 J«1,nnva 
writa( iwr, 1 19) j 

wr 1 to( 1 VP, 120) ( wxnv( 1 . J ) , 1- 1 , nnvp) 
vPlto( Iwp, 120) (wynv( 1. j ). I-I.nnvp) 
writa( IWP, 120) { wznv(1 . J ) . » - i . nnvp ) 

51 continue 

do S3 j“1,n1va 
VP 1 te( 1 VP, 119) J 

vpl te( I wp, 120) ( wx1 v(1 , J ) . l a i ,nlvr ) 
vpitel I wp, 120) ( wy i v( i , j ) , t ■ i . nivp) 
wpl te( Iwp, 120 ) (wzl v( I . J ) . 1 ai ,n 1 vp) 





42 












53 continua 

52 continue .». r<^u>4. .i: 

300 continue OF POOR OUALITV 

goto 49 ^ * 

200 continue 

K2"4 . ‘r^otaS/C (r+eta3)»»2 ♦ z3*«2) 

If (k2.eq. 1 . ) goto 47 
f •alog(4./sqrt( t.-k2)) 

e«t.+.5*(f-.5»*( 1.-k2) ♦ . ia75*(f-1.0a333)*( 1 -k 2 )»( 1 -k 2 ) 
capk«f».2S«(f*t.)*( 1.-k2) ♦ . t4*(f-1. 16666 )«( 1-h2)*( 1*k2) 
goto 48 

47 e-1. 

capkaio. * 

48 wra-fpi •dgam«2. •aqrt(r/eta3/k2 )»(capk»(2. -k2)-2. «o) 
psi"dpsln 

wz«0. 

do 250 k«t,2*kapa 

jc3"P«r * ota3*ota3 z3*z3 - 2 . •p«ota3«cos(ps1 ) 

wtercp^dpsl I •r«(r*ata3«cos(pa1 ) )/(ota3«eta3'^p«r'2. •P*eta3«cos(pa I ) ) 
wtemp*wtanip»( t . -(z3/8qrt(x3>) ) 
wz>wteRip ♦ wz 
psi"psi4'dp9t < 

250 continue 

wz»fpi ■ctg8n«wz 
If(n.oq.l) goto 210 
lf(n.oq.2) goto 230 

49 continue 
C 

C SELF-INOUCEO VELOCITIES USING SCULLY APPROXIMATION 

C 

If ( isel f .aq.O) goto 37 
if ( iwr1te.eq.2) wr 1 te( Iwr , 150) 
do 54 I«1 ,nnvr 
do 54 j ■ 1 ,nnva 

r»sqrtl xnv( i , J )»»2+ynv( 1 . J )**2) 

dwz*fpi •garat( i )«(alog{8. •r*tan1dpsin/4. )/fps2)-.25)/r 
wznv( i , J )»wznv( I , J )+dwz 

54 continue 

do 55 i"t,nlvp 

if ( 1 ,eq. t .OP. i .oq.2) ksaveaniva 
tf(1.eq.3) ksave"n1va 
do 55 J«t,kaava 

p«sqpt(xiv( i . J )»*2>ylv( I . J )**2) 

dwz'fpi •gamt I ( i )*(a1og(8. •p*tan(dpsi i/4. )/fps2)- .25 )/p 
wzlv(i,j)«w 2 lv(l,j )-*-dwz 

55 continue 

57 continue 

C 

C CONVERGENCE TEST AND NEXT APPROXIMATION 

C 

ktest'i 

i f ( ni tep . eq. I )goto 69 
savevel afacvet 
do 67 i"i,nnvp 
do 67 J • 1 ,nnva 

wxnv{ i , J )«facve1 •wxnv( i , J ) + ( 1 . -facvel ) *wxnvs( i , J ) 
wynv( I , J )«facvoi •wynv( i , J ) * ( 1 . *facvol ) ‘wynysi i . J ) 
wznv( i , J )"facvei *wrnv( I , j ) ♦' ( 1 . -facvel ) •wznvs( t , J ) 

67 continue 

CO 68 i«i,nlvp 
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o 


OO 68 J*t,n(va 

wx<v< l.i >«facv«l*wxiv( « , J ) ♦ ( 1 . 'facvel )«wxivs( i , j ) 
wy ivt t . J ) ■fscv«l «My ( v( I . J ) * ( 1 . -f acv«1 ) *wy I vs( I « J ) 
wz w( I , j )«facval *wz I v( I , j ) * ( 1 • acve I ) «wz ivs( I , j ) 

68 continue 

wconv"0. 

do 1000 i»1,nnvr 
wconv*wconv * abs(wznv( 1 ,2) ) 

1000 continue 

wconv«wconv/f 1 oat< nnvr ) 
f aevel ■savovel 
ktest«0 
do 60 1 ■ 1 , nnvp 
do 60 J • 1 . nnva 

I f (aoal wzr1v( 1 . j ) •wznva( t , J ) ) .gt. ( .02«abs(wconv) ) ) ktest*l 
' i f(kteat. eq. 1 ) goto 66 
60 continue 

do 65 Ixlinlyr 
do 65 J ■ 1 ,niva 

I r (aOa( wzl v( I , j ) «wz I vs( 1 . J ) ) .gt . ( .02*ebs(wconv) ) ) ktest* 1 
If (ktest.eq. I ) goto 66 

65 continue 

66 continue 

69 continue 

t do 203 1 ■ 1 , nnvp 

do 203 J "I, nnva 
wxnvs( 1 , J ) "wxnv( I , j ) 
wynv3( I , J )"wynv( I , J ) 
wznvs( I , j )"wznv( 1 , J ) 

203 continue 

do 204 l«l.n1vp 
do 204 J«1,nlva 
wxivsl I , j )«wx1v( 1 . J ) 
wylvs( I . j )»wylv( I . J ) 
wz1vs( 1 . j )»wz1v( 1 , j ) 

204 continue 
1f(lwplte.tt.1 )pQtupn 
wp1te( iwp, i25)ktost,facve1 
wp I tel I wp , 130) 

do 80 J "I, nnva 
wp1 to( IWP, 119) J 

wpl te( I WP , 120) ( wxnv( I , J ) , I "1 , nnvP ) 
wPl te( IWP, 120) (wynv( I , J ) , I » t ,nnvp) 
wpl te( I wp, 120) (wznv( I , J ) , I • 1 , nnvp) 

80 continue 

wpl te( twp, 140) 
do 90 J“1,nlva 
WP I te( I WP , 1 19 ) J 

wplta( I wp, 120) ( wx 1 v( I , J ) , I ■ 1 , nl vp ) 
wpl te< I wp , 120) ( wy I v( I , J ) , I • 1 ,nlvp) 

WPl tef Iwp, 120) (wz1v( I , J ) , I"! ,nlvp) 

90 continue 

return 

100 formatf/,* L0QP2: ITERATION •,l4,/,‘ GAMC, GAMT, GAMTl, GAMTIP:*) 

123 'OPmat(2{9f 10.6,/) ) 

110 fopwatl/,* POSITIONS OF RQLLEO'UP VORTICES 3f 10.6 ) 

I 19 format ( / , l4 ) 

120 format(2{9f 10.6,/) ) 

121 format!/. 14, 5x.f 10.6) 

125 format!/,* L00P2 CONVERGENCE ’.M." RELAXATION -.F5.2) 
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130 format(/.’ NEAR WAKE INDUCED VELOCITIES:*) 

MO forwatr/.* INTERMEDIATE WAKE INDUCED VELOCITIES:*) 

1*0 format!/. * SCULLY CONTRIRUTION FOR EACH TRAILER: *) 

ana 

V 








subrout Ins rollupCncall) 

Tms suoroutln* calculatos: 

m) the strengths of the trailing vorticiss in the 
near, intartaeoiate and far wakes; 
b) in addition, tha apanwise positions of the roiled up 
internodiato wake trailing vortieies, interpolat ing 
an entire new internediate wake (neat 1*1) or Just the 
first azimuthai position (ncat)"3). 

coiMson/ganuBa/gafflC( 16) ,ganit( i5),gamt 1 (3) .gant ip.eta i . eta2. eta3, 

&j 1.J2.k3,facgam 

common/parra/ iwr, ird, 1 ins, 1 ini , I im2,nl ter, iwri ta 

conmon/geom/nbidsl ,nblds.si0na.fnu,otan( iS) ,knnvr,otai(0) ,’Knivr, 

&l twist, thotadi 15),theta( i5).thotac( I4),thate0, the tOd, alphas, 
6cd0.cdk,dpsind.dpsin,dpsl id.dpsi 1 .coaff .cooff I ,c.s, blades, 
4nnvr,nnva,nnvr 1 ,nnval ,nncr,nnca,etanv( 15) .atanc( 16) , 
&nlvr.n1va.nivri .nival ,nidr,nica,otaiv(6) , Ota ic(7) , 
Sntva.ntva1.ntca.fpst.fps2 

comrnon/pos 1 t/xnv( 13, ID ) ,ynv{ 15, Ifl) , znv( 15, i8),xiv(0,5O),yiv(6.5O). 
&z1v(e,SO) 

CLEAR THE ARRAYS OF TRAILING VORTEX STRENGTHS 

1f(ncal I ,oq.3) goto 75 

do 10 l«l,nnyr 

gent( 1 )»gafflc( i+1 )»ganc( 1 ) 

continue 

do 15 i»l,3 

gamtl(1)aO, 

continue 

FIND WHERE GAMC IS MAXIMUM NEAR THE TIP (GAMT CHANGES SIGN) 

AND ROLL UP FROM THERE OUT TO TIP 

0max>O. 

do 20 l•nncr.2.•1 

lf(aos(gamc( 1 )).gt.gnax) gmax«abs(gamc( i ) ) 

1 f (gnax.eq.abs(gamc( 1 ) ) ) Jl"i 

continue 

continue 

If ( j t.gt.nnvr) wr i to( iwr .200) 

fornatt//.* Rollup: You screwed sonothing up.*) 

do 40 !■) 1 .nnvr 

gamt i (3) ■gamt i (3)'*-gamt ( 1 ) 

continue 

FIND WHERE GAMT IS MINIMUM INBOARD ANO ROLL UP ON EITHER SIDE 
OF MINIMUM INTO TWO MORE TRAILING VORTICIES 

J2«2 
J2-J 1-3 

i f (nnvr .eq. 15) j2«J 1-4 
do 60 I • j 2 . j 1 - 1 
gant i ( 2 )"gamt 1 ( 2)+gant( 1 ) 
continue 
do 70 i-a.J2-1 
ganti ( 1 )»oawt 1 ( 1 )+gamt( i ) 
continue 
gamt ip*gamt i (3) 


70 





If (ncal I .aq. 1 ) return 
C 

C MNO CENTROID OF EACH SECTION OF 8LA0E ROUEO UP 
C 

etal*0. 
eta2«0. 
ata3«0. 
do 60 

etataetat * gatntC < )«otanv( 1 ) 

80 continue 

otalaeta1/ganitt( 1 ) 
do 90 I«j2.jl-I 
eta2aota2 * gat»t{ I )«otanv( I ) 

90 continue 

eta2aeta2/gamtt (2) 
do 100 lajl.nnvr 
eta3*ota3 * garat( i )*etanv( 1 ) 

100 continue 

eta3aeta3/gaat1 (3) 

1f(nltar.gt. 1 ) return 
C 

C CALCULATE INTERPOLATION FACTORS EITHER FOR NEAR WAKE (NCALL«3) 

C OR FOR INTERMEDIATE WAKE (NITER* 1) AND INTERPOLATE NEW POSITIONS 
C 

do 110 1*i.nlvr 
1 f (otal . ge.eta iv( 1 ) ) 11*1 
1f(eta2.gs.ataiv( I )) 12«1 
If (otaS.gt.etalvC 1 ) ) *3»1 
110 continue 

terpi "(eta 1-eta lv( 1 « ) )/(etalv( 1 1+ 1 )-etaiv( 1 1 ) ) 
tarp2a(Bta2-eta1v( 12) )/(eta1v( I2-»’1 )-eta1v( 12) ) 
terp3«(eta3-etaiv( 13))/(eta1v(13>i)-ataiv(13)) 
do 140 J ■ 1 ,ni va 

X 1 v1 “X 1 v( 1 1 . J ) torp1*( X 1v( 1 14-1 , J )-x 1 v( 1 1 , J ) ) 

xiv2"xlv( 12. J ) terp2*(x1v( 12+1 , J )-x1v( ! 2. J ) ) 
xlv3axiv( 13. J ) ♦ tarp3*(x1v( lO+l , J )-xlv( 13, } ) ) 
ylvl ay 1 v( 1 1 . J ) ♦ terpi «(y 1 v( 1 1 , J ) -y 1 v{ 1 1 , J ) ) 
y lv2»y lv( 12. J ) ♦ torp2»(y 1 v( 1 2* 1 , J )-y I v( 12, j ) ) 
y1v3«y 1v{ 13. J ) torp3« { y 1 v( 1 3> 1 , J ) -y I v( 1 3 , J ) ) 

z1viailv{ 1 1 . J ) terpi*(zlv( 1 1+-1 , J )-z1v( 1 1 . J ) ) 

z1v2"Z 1 v{ 12. J ) ♦ terp2»(z1v( 12* 1 , J ) -z1 v( 12 . J ) ) 

2lv3asiv( 13, J ) + terp3*(z1v(13-*-1,J )-z1v(13,j)) 

X lv( 1 . J )«xlvi 
xlv( 2, J )*xiv2 
x1v( 3 , J )«xlv3 
y lv( 1 , J )«y ivl 
ylv(2, J )ay1v2 
y 1 v( 3 , J )ay 1v3 
Zlv( 1 , J )»ziv1 
z 1v( 2 , J )"zl v2 
ziv(3, J ).aziv3 
140 continue 
nlvr«3 
nlcr«nivr*l 
nlvrl«ntvr-l 
return 

73 continue 
x«0. 
y»0. 

Z»0. 
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/ 




155 


156 


157 








do 155 1-2.J2-1 

xaxnv( 1 .nnva)*genit( 1 )♦« 
yaynv( I .nnva)«gent( I )*y 
x«xnv( I ,nnva)*gant( 1 )+z 

continue 

Klv( 1. 1)>x/gamt1( 1) 
y1v( 1,1 )ay/cantt( 1 ) 
zlv( 1, 1)az/gamt1( l) 
xao. 
y»0. 

*■ 0 . 

do 156 laJ2,J1-1 
xaxnv( 1 ,nnva)*g&mt( 1 )*x 
yaynv( I ,nnva)agarat( 1 )+y 
raxnv( 1 ,nnva)*gamt( 1 )-*-x 
continue 

x1v(2, t)ax/ganit1(2) 
y1v(2, 1 )ay/gamti(2) 
Ziv(2, 1 )az/gant1 (2) 
x«0, 
y»0. 

X"0. 

do 157 lajt.nnvp 
xaxnv( 1 ,nnva) agarntt t )*x 
yaynv( 1 .nnva)«gamt( I )+y 
xaxnv( 1 ,nnva)*gamt( 1 )+x 
continue 

x1 v( 3 , 1 ) ax/gamt Ip 
y1v(3, 1 )ay/gamt1p 
X 1 v< 3 , 1 ) az/gamt Ip 
return 
end 





c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


300 


c 

c 


11 


to 

c 

c 

c 


subroutine vlnd 
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This routine calculates the InOucod velocities on the 
blade (VINOS), the near wake (VlNON) and the Intermediate 
wake (VINOI), due to a sagnent element of endpoint coordinates 
xl.y1.zi and x2.y2.z2 and of strength gm. VINOS only calculates 
the tnfluonco coefficients cw( 1 , I ,n) , which give the velocity 
induced at spanwlse position I by a trailing vortex of strength gm«1 
located at spanwlse position I (sumn«d over all azimuthal positions) 
in waxe n (noar«l, lntorn)edlate>2. fara3). These are used in LOOPl. 

VELOCITIES INDUCED ON 8LA0E 


coauson/parn/iwr. Ird, I Ims, I Ini . I lra2.nl ter, Iwri te 

common/geom/nolds 1 . nb I do. a 1 oma. fmu.et an( iS ) .knnvr.eta 1 (6) , kni vr . 
ftl twist. thetadf IS ) . thota( IS) . thotac( 14). thotaO. the tOd. alphas, 
tcdO.cdk .ops lnd.dpsln.dps lid, dpsi 1 ,coof f ,coef f 1 .c. s.blades, 
&nnvr,nnva.nnvrl .nnval ,nncr.nnca.otanv( iS) ,etanc( 16). 

&nlvr .niva.nfvri .nival ,nier,n1ca,ataiv(C),etaic(7) . 
ftntva.ntval .ntca. fpsi ,fps2 

common/vindo/lfar.nfar.xt ,x2.yl ,y2.zl .z2.x,y.z.facvel .gra.ux.uy.uz 
C 0 mnon/p 08 it/xnv( 15, 10) ,ynv( 13, 18) ,znv( 13, 18) ,x1v(6.S0) ,y 1v(6,50) . 
Az1v(6,50) 

common/valoc/wxnv( 13,18) ,wynv( IS. 18).wznv( IS, 18),wx1v(6.50). 

&wy1v( 6.S0) . wz 1v(6, SO) .facgoom 
COfflmon/coef /cwx( IS. IS, 3) ,cwy( IS. 15,3) ,cwz( 15,13,3) 
comraon/garama/gancC 16) ,gant( IS) .gant 1 (3) .gant lp.eta(3) . 

&J 1 .kk(2) , faegara 
entry vlnab(nl .n2. 1 ,n) 

If(oni.ne.l) wrl te( lwr.300) 

formate/,* VINOS: What are you doing? gn Is not equal to one.") 

tsk 1p"0 

do 10 1*<,nnvr»1 
x«(etanv( I )*etanv( l+t ))/2. 
y»0. 
z-0. 
cc« 1 . 
as>0. 
n3»1 

If (n1 .eq.n3.or.n2.eq.n3) iskfp*1 
do 10 j«i,nbids 
If ( iskip.eq. 1) goto 11 
call wxyz(fps2) 

cwx( 1 . 1 .h J«cwx( 1 , 1 ,n) •*> (ux^cc uy*ss) 

cwy{ 1 , 1 ,n)»cwy( 1 . I ,n) ♦ (uy*cc • ux*ss) 

cwz( 1 . 1 ,n)«cwzO . 1 .n) ♦ uz 

continue 

1 sklp*0 

sav*x 

x«x»c - y«s 
y«y»c ♦ sav"s 
sav"cc 

cc»cc«c - ss*s 
ss"ss*c ♦ sav*s 
continue 
return 


VELOCITIES INDUCED ON NEAR WAKE 


entry v1nOn(n1 , n2 . loop) 


o o o 
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lsk<p*0 i 

r.f»nnva ! 

<f(nl.eq.O) nf«nnva-t | 

if I loop.eq. 1 ) nf«1 j 

do 20 i • t , nnvr \ 

do 20 J"t ,nf I 

x"xnv(t,J) i 

y«ynv(l.J) | 

X"znv(l,J) i 

CC«1. j 

ss>0. j 

■ n3"( J - 1 >*nnvr ■*>1 

tf(n1 .«q.n3.or.n2.aq.n3) taktp*l I 

do 20 k<> 1 .nblds ! 

( f ( 1 ak Ip.eq. 1 ) goto 21 
cat I wxyz( fpa2) 

wxnv( 1 . J )»wxnv( 1 , j ) ♦ (ux«cc ♦ uy«sa) 

• wynv( 1 , J >"wynv( 1 . J ) ♦ (uy«cc • ux*35) 
wznv( I , J )«wznv( 1 , j ) ♦ uz 
21 continua 
tsklp-O 
sav«x 

x»x»c * y*9 
yay»c ♦ aav*8 
aav^cc 

cc»cc*c * 3a«9 

93"39>C ♦ sav*9 

_ 20 continua 

/ I f < loop.eq. 1 ) roturn 

If(nl.no.O) return 
taktpaO 

do 26 1«1,nnvr 
x*xnv( 1 .nnva) 
y»ynv( 1 ,nnva) 
z»znv( 1 .nnva) 
cc* 1 . 

99-0. 

r-9qrt(x*x ♦ y*y) 
ri-sqrt(x1-x1 y1-y1) 
if (abs(r*p1 ) . I t. .025) I9k1p-1 
do 26 k-1.nb1d9 
1f( isktp.aq. t) goto 26 
call wxyz(fps2) 

wxnv( 1 .nnva)-wxnv( I .nnva) ♦ (ux-cc * uy-99) 
wynv( 1 , nnva ) -wynv{ t .nnva) ♦ (uy-cc * ux»sa) 
wznv{ 1 .nnva )-wznv( 1 .nnva) uz 
28 continue 

lak ip«0 
aav-x 

x-x-c - y-3 
y-y-c ♦ aav-9 
sav-cc 

CC-CC-C • 9B-9 
33-39-C t 3av-9 
26 continue 

return 

VELOCITIES INDUCED ON INTERMEDIATE WAKE 
entry vindl(nl,n2) 







•3‘ r’ 

rr- If 


ORIGINAL ® 

OF POOR QUALTTf 



Isk1p«0 
nf • I 

1f(n3.oq.O) nf»a 
do 30 I't.nlvr 

If(l.flq.t) ksavo«n1va 
1f(<.eq.2) ksava«n1va 
<f(f.eq.3) ksavsanlva 
do 30 j*nf,kS8vs 

K«X(V( 1 . J ) 

y»yW( i.i) 

*•! lv( ♦ , j J 
cc« 1 . 
sa*0. 

n3»{ J - 1 )»nlvp nnva*nnvr i 
tf (n1 .eq.n3.or.n2.aq.n3) 1skip«1 
do 30 k>1,n01d3 
4f (isktp.eq. 1 ) goto 31 
call wxyz(fps2) 

wxlv( 1 , J )»wxlv( 1 . J ) ♦ (ux*cc ♦ uy«ss) 
wytv( t . j )«wylv( 1 . J ) + (uy*cc - ux«sa) 
wzlv( I . j ) ■wr w( 1 , j ) ♦ uz 
31 continue 
lsklp-0 
sav"x 

x"x»c - y»9 
y«y*c ♦ sav*9 
aav'cc 

CC«CC«C • S3*3 

ss«ss«c sav*s 
30 continue 

lf(n2.no.O) return 
lsklp»0 

do 33 lal.nlvr 
x»xlv( 1,1) 
y*y lv( 1,1) 
z«z lv( 1,1' 

CC" I . 

ss-0. 

r«aqrt(x*x y«y) 
r2»sqrt(x2»x2 ♦ y2*y2) 
lf(aDS(r-r2). It. .025) laklp-l 
do 35 k*i.n0ld3 
If ( l3klp.eq. 1 ) goto 36 
call wxyz(fps2) 

wx I v( 1 , 1 ) "wx I v( I , 1 ) ♦ ux»cc ♦ uy*S3 
wy lv( 1 , 1 )»wy lv( 1 , 1 ) ★ uy»cc • ux»cc 
wzlv( I , 1 )«wzlv( I . 1 ) ♦ uz 
36 continue 

I sk lp*0 
sav>x 

x«x«c - y»s 
yay«c ♦ SaV*S 
sav-cc 

CC"CC*C • 3S*S 
SS«33»C ♦ sav*s 
33 continue 

return 
end 







tubrouttna wxyx(fpa3) 

PO'nt X.y.,. It I, only C.liw 5Mr^?' cm on «ny 

rail. 4 10 

yct».5«(y1*y2) 

iCt«.S*(2Ux2) 

aXKax-KCt 

Oyyay.yCt 

a**«x-zct 

fia»»qpt((x1-x)*.aa,y, 

Osx».5«(x2-x1 i •'^-*3. lo. . l*-3) goto 100 

a«y».3»(ya-y t ) 

Os**.5«(z2-xi) 

OS2-09X..2 ♦Osya ♦dax-.a 
fva*fos2*fp32*ds2 

rmax2«400. *032 

r02*dxx«*a ♦ dyy**2 ♦ dzz**2 

r03«r02*sqpt(r02) 

dSAxadSZ«dyy . dSy*dZZ 
C9myd3x*dzx - dsx-dxx 
Osiiiz«dsy*oxx - dsx»dyy 

o»'"y*a ♦ d5mz..2 
jr(dsm2. l«, . ie.)3) goto 100 
fvds-fv«/dsn2 

»0«l. 

1Mr02.gt.pmax2) goto 120 

«lpaaaatpna2 - a«a 
< f (abs(atpaa ) . la. . ta- 15 ) ootn 
«qta-3qrt(aoi( 1. ♦ 2. *a ♦°a?Dha3?i 
*qaa.sqrt(ab3(i. - 2..I I f ! 

|f(SQla. 1 t.. 001 ) 9qla..00l 
^ * 0 *®- • t. . 001 J sqiao.oot 

^ (..p..2-.,/.daa 

ooto 120 
cont 1 nue 

*0" >./(( 1 . *a*a)««2) 

fact»* 10 *fpl*grn« 2 . /p03 

fact«fact/( l.+fvds) 

ux»fact*d3mx 

uyafact«asmy 

uz"fact«dsfflz 

Paturn 

cont Inua 

ux»o. 

uy*0. 

U2"0. 

paturn 

and 
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subrouttno (ntgr(ktest) 
c 

C This routlno Integrates the wake geometry using the InOucad 
C velocities returned by LOOPS. Tangential and radial velocities 

C are not used; the wake is integrated directly from WKnv, wynv 

C and wsnv. The velocity of a wake nods In travelling from 

C position l,j to Itj-f-i IS assumed to be tne average of the 

C velocities at t.J and l.Jti. 

C 

comnon/parm/lwr. ird, I iffls. I Imt . I lm2.nl ter, Iwri te 
comnon/geom/nbldsi .nulds.signa, fmu.otani 15) .knnvr.etal (6) .knlvr, 

&l twist . thetadC 15) . theta ( 1S) . thetacl 14 ) , tnetaO, t ha tOd, alphas, 
AcdO.Cdk .dps Ind.dpsln.dpsl Id, dps 1 1 .coef f .coaf f 1 ,c. s. blades, 
&nnvr,nnva,nnvr1 .nnval .nner .nnca.etanvC 15) ,etanc( 16) , 
&n1vr,nlva.nlvri,nlvBi.nlcr,ntca,eta1v(6),etaic(7), 

&ntva,ntva1 ,ntca. fpst , fps2 

common/poslt/xnv( IS. 18hynv( IS, IS) ,znv( IS, 16) ,x1v(0,SO) ; 
&y1v(6.50).2lv(6,50) 

eomnon/veloc/wxnv( 1S. I8).wynv( IS, 18),wznv( IS, 16) ,wxlv(6,50) , 
&wy1v(6,50) . wz1v(6,50) . faegaora 

COmmon/gamma/gamcC 16) ,gcm't( IS ) .garet 1(2) .garat Ip.otaf 2 ) ,k 1 ,k2 ,k2. faegam 
common/ I ntgrsave/xnvs( 15,18) ,ynvs( IS, 18 ) ,znva( 1$, 18 ) , x I vs(6 , 50} , 
&ylvs(6.50).z1vs(6,50) 
common/speclal/ I tgr 

ktest "1 

If (nl tor.eq. 1 ) goto 202 
do 10 l• 1 ,nnvr 
CO to j«l,nnva 
xnv(l.j)»0. 
ynv( 1 , j )"0. 
znv(1,J)»0. 

10 continue , 

do 20 I ■ 1 , nnvr 
xnv( 1 , 1 )«otanv( I ) 
ynv( 1 . 1 )»0. 
znv( 1 , 1 )«0. 

ph1«0. 

do 20 J"1 ,nnva-1 
phi 1 *00 1 tops In 

wr»('-xnv< 1 , J ) * 008 (phi )*wxnv( 1 , J*1 ) •cos(ph1 1 )»wynv( 1 , J )*sln(phl )♦ 
6wynv( I , J-M )«sln(phl 1 ))/2. 

ut«( -wxnv( I,j)*s1n(ph1) «wxnv( I , )♦ 1 ) •s1n(ph1 1 )*wynv( I , J )*cos(ph1 )♦ 
6wynv( I , J+1 )*cos(phl 1 ) )/2. 
r«sqrt(xnv( I , J )»«2 + ynv(l,J)»*2) 
r1«r»apsin»wr 

phi t»phl+dp3ln*(2. •wt/(r+r1 ) )>dps1n 
xnv( I . Jt-1 )»ri - 008 ( 00 1 I ) 
ynv( 1 ,J^1 )"r 1 - 5 ln(phl 1 ) 

znv( t , )“znv( 1 , J ) ♦ dpsln-(w 2 nv( 1 ,J )*wznv( I , J+1 ) )/2. 
phi -phi 1 
20 continue 

oo 20 l•l,nlvr 
oo 30 J • I .niva 
xlv(1,j)-0. 
yiv( I . J )»0. 
ziv( 1 . J )«0. 

20 continue 

nca 1 1 “3 

call rol lup(ncal 1 ) 

1f{ itgr.ne. 1 .or.nnva.gt.2) goto 1033 
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ao 1030 l•k1,nnvr 
xnv( 1 ,nnva )»xiv(nlvi*.1) 
ynv( 1 ,nnva)«y1v(n1vp, 1 ) 
znv ( 1 , nnva )»z1v(n1vp,t) 
continue 
continue 
zIvCnIvr, 1 )».025 
do 40 l•l,n^vr 
phi •atan2(y I v( l.l).xlv(l.l)) 
do 40 J • I , nlwa- 1 
phi l«phl^dpBl I 

wp*( wx lv( 1 , J )*cos(phl )+wxlv( I , J+ t )*co5(ph1 1 )*wy1v( I , J )»sin(phl )♦ 
ftwylv( I , J-M )«sln(phl t ) )/2. 

wt«( -uxl v( I . J )»s1n(ph1 )-wxlv( l,J*1)»9ln(phl1 )*wy ivl 1 , J ) *cos(ph1 )* 
&wylv( I , J )*cos(phl t i )/2. 
r»9qpt(xlv( I ,J )»»2 ylv(l,J)*«2) 

r1»p+dp9l l«wp 

phi l-phl+dpsi I♦dps1 l•{2.•Mt/(r♦^1 ) ) 
xlvll , J^i )»ri*cos(pnl I ) 
ylv( I , j*l )«r1»sln(ph1 1 ) 

zl v( I . )«z1v( I , J ) ♦ Opal I*(wz1v( 1 , J )^wz1v( 1 , J*l ) )/2. 

phl«phl I 

continue 

savegaom* r acQaon 

do 41 l"1,nnvr 

do 41 J • 1 , nnva 

xnv( I , J )«xnv( I , J )»facg6om xnvs( I . j )•( 1 . -facgeoni) 

ynv( I , J )»ynv( I , J )*facoaom ynvs( i . J ) •( 1 . -facgoom) 

znv( 1 , J ) «znv( I , J ) *f Bcgoom ♦ rnvs( 1 . J )•( 1 • *^acgoom) 
continue 

f acgeoniBsavegeom 
do 42 1»l,nlvr 
do 42 Jal.niva 

xlv( I , J )»xlv( I , J )*facgeom ♦ xivs( I , J ) *( 1 . *facgoom) 
y1v( I , j )«y1v( I ,J )*facgeoni ♦ y lva( 1 . j ) •( l . -facgoom) 
ziv( I , j )»z1v( 1 , J )*facgeom * zlva( 1 , J }•( I .*facgaam) 
continue 
facgeoniaaavogeoni 
Ktaat*0 

do 43 jBl.nnva 
zconv"0. 

do 1000 ll■l,^nvr"1 
zconvazconvi-abaCznvC 1 1 ,J )) 
continue 
zconv»zconv/nnvr 
do 43 I ■ I ,nnvr- I 

I f (aba (znv ( I . J )»znva( I . J ) ) .gt. ( .02»zconv) ) kteat> I 

If (ktcat .eq. I ) goto S3 

cent Inue 

do 44 i»i.nlvr 

do 44 j»i ,nlva 

I f (aba( z I v( l,J)*z1v3(1,J)).gt. -02*aba(z I v( I.J))) ktost»l 

I f (k test . eq . I ) goto 53 

continue 

continue 

cont inue 

do 202 1 ■ I , nnvp 

do 202 J “ 1 . nnva 

xnva( 1 , J ) •xnv{ l , ] ) 

ynva( 1 , j )»ynv( 1 . j ) 


♦ xivs( I , J )*( 1 . -faegeom) 
•k ylva( 1 ,j)*( 1. -facgoom) 

♦ zlvs( 1 . J )•{ I .-faegaom) 


•Ot* ( .02»zconv) ) ktaat»1 


,gt. -02*aba(z1v( I . j ) ) ) ktost»l 








znv8< i . J )"znv( I , J ) 

202 contlnua 

do 204 (Bl.nlvr 

f f ( 4 .«q. f .or. I .oq.2) l(sava"nlv« 
tf(l.«q.3) k8ave*n1va 
do 204 Jal.kaav* 
k4vs( 4.J)*x1v( 4.J) 
y4va( 4 . j )»y tv( 4 ,J ) 

Z4vs( 1 ,j )«z4v( I ,J ) 

204 cont4nuo 

4f ( 4wr 4 to. 1 1 . 4 ) goto SO 
wr I ta( 4wr. 400) ktost. facgeom 
do 4S J"4,nnva 
wr4 to( Iwr , 409 ) J 

wr 4 to( 4wr , 4 40) ( xnv( 4 , J ) , 4*4 ,nnvr ) 
wr4 to( 4wr , 4 40) (ynv{ 4 , j ) , 4*4 ,nnvr ) 
wr 4 to( 4wr , 4 40) (znv( 4 , j ) , 4 «4 ,nnvr ) 

43 cont4rHjo 

do 46 J"1.n4va 
wr 4 to( twr , 409) J 

wr4ta(4wr.440) (x4v(4,J),l»4,n4vr) 
wr4to( 4wr, 110) (y 4v( 4 , J ) , 4« 1 ,n4vr ) 
wr 4 te( Iwr, 1 10) ( z 4 v( 4 . J ) . 4 • 1 ,n4 vr ) 

46 continuo 

50 return 

400 fornaK//.* INTGR CONVERGENCE •,44.* RELAXATION ".F10.S) 
109 foPmat(/. 44) 

4 40 fOPiTiat(2(9f 10.6./)) 
end 
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subrout I n« loopMktast) 

This routina calculatas a new distrfbutton of bound circulation 
gamc(l) consistent with tho current wake geomotry as given by 
cw(i.l.n). by Iteration. 

INFLUENCE COEFFICIENTS CW( 1,1.1) 

common/veloc/wxnv( 15, 18),wynv( 15. 18),wznv( 15. 18) ,wx1v(6.50) . 

&wy1v(6.S0) ,wz1v(6.50) .facgeom 
comiDon/parm/ Iwr . ird, 1 1ms, 1 1ml , 1 lni2, niter, iwrlte 
conmon/geoRi/nblds 1 .nblds.slgma, fniu.etan( IS) .knnvr .otal (6 1 .knivr , 
ft I twist , thotad( 15) . theta( 15) , thetac( 14 ) , thetaO, thetOd. alphas. 
ftcd0.cdk.dpslnd.dps1ri.cps1 Id.dpsi 1 .coef f .cooff 1 .c, s. blades, 
ftnnvr.nnva.nnvrl .nnval ,nncr.nnca.etanv{ tS) ,etanc( 16 ) , 
ftnl vr .niva.nivr 1 ,n1 va 1 .nicr .nica.eta I v(6) , stale! 7 ) , 
ftntva.ntval .ntca, fpsi . fps 2 

conunon/gamma/ganci 16 ) .gamt ( 15) .garnt 1 (3 ) .gant ip.etal ,eta2 .eta3. 
ftj 1 ,kk(2 ) . faegam 

common/ V I ndo/ 1 f ar.nfar ,x 1 , x2.y 1 .y2. z 1 ,z2. x.y ,z, facvel .gra.ux.uy.uz 
conmon/posi t/xnv( 15. 18).ynv( 15. 18), znv( 15, 18 ) , x1v(6,S0) . 
ftylv(6.50) .Zlv(6.50) 

common/coef/cwx( 15. 15,3) ,cwy( 15. 15,3) ,cwz( 15,15,3) 
common/sal f /I set f 

comreon/extraspec/fps3 
dimension games! 16) 

common/ 1oopsave/wxv( 16) ,wyv( 16) ,wzv( 16) 
data pi .twopi . fp 1/3. 14159,6.28318. .079577/ 
feal k2 
do 5 1 ■ t . nnvr 
do 5 J ■ 1 , nnvr 
do 5 n« 1,3 
cwx( 1 , J ,n)"0. 
cwy( I , J ,n)«0. 
cwz( I , J ,n)aO. 

5 continue 
dps2"fps2 
do 10 1*l,nnvr 
fps2*fps1 

If ( 1 .eq.nnvr-1 ) fps2"dps2 
do 11 J»l,nnva“1 
x1»xnv( 1 , J ) 
y1»ynv(1.J) 
z1"znv( 1 , J ) 
x2«xnv( 1 . j* 1 ) 
y2»ynv( 1 . J-M ) 
z2"znv( 1 , )•*■ 1 ) 
ni«( J - 1 ) «nnvr ♦ 1 
n2«nl ♦nnvr 
gm« 1 . 

call vlnab(n1,n2, 1.1) 

11 continue 
10 continue 
fps2«cps2 
pn ib- twap1/b1 ades 
CO 20 1»1,nlvr 

If(l.eq.l) ksave>n1va 
lf(i.oq.2) ksavo«nlva 
1f(i.eq.3) ksavs“n1va 
do 21 J»l.kaavo-1 
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xi"xiv(r,j) i 

y1»y«v(1,J) ; 

z1«z1v(i.J) ; 

x2*xW( t . J + 1 ) 

y2«yW(1,J«'1> j 

z2*z1v( « , J*1 ) i 

ni«( J - I ) •nl vr ♦ nnvr*nnva +1 | 

n2*nl+’nivr I 

aps2«fps2 j 

pnl l»f loat(nnva- 1 )»Ppain ♦ f loat( J • 1 )*PpsH 

If (pm 1 .gt.phlb-dpstn) fp92afps3 • 

gm» 1 . 

call vlndb(n1 ,n2, I .2) 
fps2*dps2 
21 continue 

20 continue 

n«3 

If(lfar.eq.l) goto 401 • 

dz"Zlv(n1vp.niva)*ziv(nlvp,n1va«1 ) 

x1«x1v(nivp,n1va) 

yi“y Iv(mvp.nlva) 

z I “z 1 v(n I VP , n1 va ) 

ps1«atan2(y1 ,x1) 

p«aqpt(x1«x1 ♦ y1*yl) 

do 30 J» 1 ,nfap 

ps I "ps l-i’dps 1 1 

x2*p*cos(psl ) ' : 

y2»p*sin(ps1 ) 

z2»z1+dz 

nl“J ♦ nnvB»nnvp nlva*nlvp 
n2«nl+1 

gm«l. ; 

call v1ndb(n1 ,n2, 1 .3) j 

x1«x2 I 

y 1 -y 2 I 

z1»z2 

30 continue ! 

goto 301 I 

401 continue * 

f kapa-phib/dps 1 1 
kapa" 1 f 1 x( f kapa * .5) 
do 505 kl«1.nlvp 
Xl«x1v(k1 .ntva) 
yl«ylv(k1 ,ntva) 

z 1 «z W(kl ,nlva)-*-(zW(k1 ,n1va)-ziv(k1 .nlva-kapa)) 
dgam" 1 ,/(zlv(ki ,niva)-zlv(kl , ni va*kapa ) ) 
p«(sqpt( X 1 *x 1+y 1 *y 1 ) ) 
do 201 l-l.nnvp-l 
x»(etanv( 1 )+etanv( 1+1 ) )/2. 
y»0. 
z«0. 

z3"ab9(z1-z) 
ota3»sqpt(x*x+y»y ) 
k2«4, •r*ota3/( (r+ota3)*»2+z3»*2) 
if (k2.aq, 1 . ) goto 202 
f «aIog(4 ./aqpt ( 1 . -k2 ) ) 

e-1.^.5»(f-.5)«( t.-k2) ♦ . 1fl75*(f-1 .08333)*( 1 . -k2)»( 1 . -k2) 
capk-f ♦.25*( f -1 . )»(1.-k2)+.14»(f-1. 1 3666)*( 1 . -k2)*( 1 . -k2 ) 
goto 203 
202 e- 1 . 
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eapk«iO. 

continue 

wr«-fp1 •Ogam*2. •sqrt(r/eta3/k2)*(capk*(2. -k2)-2. 

ps I ■Opt In 

wz«0. 

do 204 k«1,2*kapa 

x3»r«p ♦ eta3*eta3 ♦ z3*z3 - 2. •r*eta3*cos(psl ) 

wteiflp»dpil 1 •p»(r*eta3*coB(P9l ) )/(ata3*ota3*p*r-2. •r««t83*co»(p8i ) ) 

wtcff>p"wtomp*( 1 . '(z3/sdrt( m3 ) ) ) 

wz«wz<»wtefflp 

psi«pst-»dp8l 1 

continue 

wz"wz * fp I •dgam 

cwm( 1 ,ki ,n)«cwx( I ,kf ,n)+wr 

cwz( 1 ,k 1 ,n) ■cwz( 1 ,k 1 ,n)>wz 

eont tnue 

continue 

Loop on gafflc(l). First clear array of btade induced velocities. 
Then calculate new wxv,wyv,wzv from Influence coefficients. 

Then new game. Back to beginning until converged. 


continue 
I f ( Iwr 1 to. oq. 2) 
ktefnp»0 

do 200 k* 1 , I Imt 
savegamaf Began 
cont inue 
do 35 ••l.nnwr 
WKV( 1 )«0. 
wyv( I )»0. 
wzv( 1 )«0. 
cont Inue 
do 36 J«1,nnvr 
do 37 l■ 1 .nnvr- 
wxv( 1 )"wxv( I ) + 
wyv( I )"wyv( 1 ) ♦ 
wzv( I )«wzv( 1 ) ♦ 
cont I nue 

If (ktcnp.eq 
cont tnue 
do 40 J»l,nivr 
do 41 i»1,nnvr- 
wxvl 1 )<»wxv( I ) ♦ 
wyvf 1 )«wyv( I ) ♦ 
wzv( 1 )»wzv( I ) 
cont inue 

If (ktenp.eq 
continue 
do 50 J*1.nlvr 
oo 51 l•1,nnvr- 
wxv( 1 ) «wxv( 1 ) ♦ 
wyv( I ) «wyv( 1 ) ♦ 
wzv( I )«wzv( 1 ) 
continue 

I f ( ktomp.eq 
continue 

1 f (ktoinp.oq 
1 r ( I nel f . eq.O) 
do 56 !• 1 ,nnvr 


wri te( Iwr, 304) 


cwx( 1 , j, 1)*gant(J ) 

Cwy( 1 ,J , 1)*oamt( J ) 

Cwz( 1 , J , 1 ;*garat( J ) 

. 1) wr1te( iwr,3i0) (wzv(k) , k« 1 .nnvr ) 


cwx( 1 , j,2)*oarat1( j ) 
cwy{ 1.J,2)«gamtl(J ) 
cwz( I ,J .2)*garat1(J ) 

.1) wr1te( Iwr .310) (wzv(k) .k« 1 .nnvr) 


cwx( I .) ,3)»gamt1( J ) 

Cwy( I , J .3}»gamtl(j ) 
cwz( I , J . 3)»gamt 1 ( j ) 

,1) wrlte( Iwr.3l0) (wzv(k) ,k*»1 .nnvr ) 

. 1 J goto tool 
goto 35 
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70 

C 

C 

c 


80 

81 


90 

94 

c 

c 

1001 


52 


91 


200 

205 


311 

310 

304 
300 

305 
400 


r*«tanv( l ) 

dwz'fpi «gant( 1 )«(alog(a. •r*t«n(dp*ln/4. )/fpo2)- .25)/r 
wzv( 1 )»wzv( 1 (♦owz 
continue 
continue 
gamc( l )«0. 
gamc( nncp)»0. 

00 70 1*1 ,nnvr-1 
wz«(wzv( 1 

wy«wyv(i) ♦ (otanv( 1 )^otanv( 1+1 ) )/2, 

t lam«wz/wy 

f I amda«atan( t Ian) 

u*9qrt(wz*«2 ♦ wy*«2) 

atp'flanda - tnetac( 1 ) 

if (alp.gt.elphail alp«alpna« 

game( M>i )ap 1 *p 1 * 8 igm 8 *u«alp/blados 

continue 

gamc(nnvp)"gafnc(nnvr)*cobf f 1 

CONVEROeNCE TEST 

kta»t*1 

If (k.ag. 1 } goto 94 

ktest*0 

do 80 l»i.nncr 

1 f (et 3 s(gafflc( 1 )-gaincs( 1 ) ) .gf. .01*ab5(0Bi*c( 1 ) >) ktost-1 
1 f ( k test . eq. 1 ) goto 81 

continue 

continue • 

do 90 1 ■ 1 , nncr 

gamc( 1 )«ganc( 1 )»facgara ♦ gamca( 1 )•( 1 . -facgam) 

continue 

continue 

If (ktest.eq.O) kternp>1 
if (ktest.eq.O) goto tOOO 
continue 

If ( 1wrlte.ne.2.and.ktest.ne.0) goto 52 
If ( Iwr 1 te. 1 1 . 1 ) goto 52 
wr I tp( Iwr ,3 1 l ) k 

wr 1 tel Iwr. 3 10) (garac( 1 ) , 1»1 ,nnvr) 

wr I te( twr . 3 10) (wxv( 1 ) , I • 1 ,nnvr ) 

wr 1 to( Iwr . 3 10) (wyvl I ) , 1" 1 .nnvr ) 

write! iwr, 310) ( wzv( 1 ), i ■ 1 , nnvr ) 

continue 

oo 91 l«1,nncr 

gamcs( 1 )*gafflc( 1 ) 

continue 

ca I I rol 1up( 1 ) 

If (ktest.eq.O) goto 205 

continue 

continue 

facgamasawogam 

I f (ktest .oq. 1 )wrl te( Iwr . 300) 

If (ktest.eq.O. and. Iwrlto.gt.O) wri te( Iwr ,305) k 

format! / . 14 ) 

format (2(9f 10.6,/) ) 

format!//. ‘ L00P1 RESULTS"./) 

format!//,* L00P1: NO CONVERGENCE") 

format!//. ■ L00P1: CONVERGED IN ". 14 ITERATIONS" ) 

return 


■ > ) 

SAMPLE OUTPUT 


••• INPUT ••• CASE # : 1 

NUMBER OF BLADES 2 
SIGMA 0.04640 

MU 0« 00000 

PITCH ANGLE OISTRIBUTIONlOO) : I7.SOOOO tS. 30000 14.20000 12.85000 10.90000 

9.25000 8.70000 8.15000 7.93000 7.71000 7.4SOOO 

7.27000 7.0S000 

stall angle 0.20000 

CO- . 0140*0. 500»ALPHA*«2 

MAX. NUMBER OF ITERATIONS FOR L00P1 AND L00P2: 20 200 

CORE SIZE FOR NEAR WAKE INBOARD: 0.050 TIP: O.OtO 

NEAR WAKE DEFINITION: ( 13. 8). 0.10000 0.25000 0.35000 0.50000 

O.6SCO0 o.aoooo O.BSOOO 0.90000 0.92000 0.94000 

0.960CO 0.98000 1.00000 

INT. WAKE DEFINITION;! 6.3S). 0.10000 0.25000 0.60000 0.85000 

0.95000 1.00000 • 

CORE BURST TO 0.05 
TIP COEFFICIENT 0.00 


in 

VO 






> 




■> 


Mtlii RESULTS -- CONVERGENCE: O 


CI» -0.004533 CP- 0.000337 
GAHC.WXC.WYC.WZC: 


-0.00628 -0.01043 -0.01269 -0.01391 -0.01533 -0.0172S -0.01989 -0.02106 -0.02119 -0.01932 -0.01443 0.00000 CAMC 

0.00125 0.00365 0.00491 -0.01000 -0.02330 -0.05013 -0.05587 -0.04742 -0.04134 -0.03B03 -0.03172 -0.02029 UXC 

-0.00037 -0.00123 -0.00233 -0.00287 -0.00276 -0.00262 -0.00225 -0.00140 -0.00080 -0.00034 -0.00005 0.00085 WVC 

0.02283 0.03163 0.04364 O.OS6B4 0.06046 0.05373 0.04163 0.03538 0.03439 0.04172 0.06219 -0.05524 WZC 


WAKE GEOMETRY: 


1 


0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

O.OCO 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0. 10000 

0.25000 

0.35000 

0.50000 

0.65000 

o.aoooo 

0.85000 

0.30000 

0.92000 

0.94000 

0.96000 

0.9SOOO 

1.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

O.OOCOO 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

2 

9.803 

9.889 

0.891 

9.913 

9.941 

9.958 

9.960 

9.977 

9.982 

9.991 

10.007 

10.006 

10.003 

0.03999 

0.2SC54 

0.35094 

0.49951 

0.64669 

0.79266 

0.84057 

O.B912S 

0.91248 

0.93420 

0.05S&S 

0.97439 

0.99061 

0.00398 

0.00661 

0.00868 

0.01123 

0.01212 

0.01194 

0.01158 

0.00949 

0.00384 

0.01046 

0.01271 

0.00098 

-0.00981 

3 

19.464 

19.739 

19.755 

19.812 

19.872 

19.911 

19.911 

19.949 

19.961 

19.983 

20.035 

20.011 

19.997 

0.09991 

0.25101 

0.35182 

0.49899 

0.64341 

0.76559 

0.83145 

0.86280 

0.90572 

0.93072 

0.9Q796 

0.96817 

0.97260 

0.00830 

0.01438 

0.01865 

0.02361 

0.02537 

0.02513 

0.02550 

0.02170 

0.02260 

0.02373 

0.02143 

0.00302 

-0.01650 

4 

29.273 

29.633 

29.663 

29.731 

29.614 

29.871 

29.867 

29.927 

29.948 

29.080 

30.076 

30.017 

29.969 

0.09332 

0.25149 

0.35270 

0.49851 

0.64019 

0.77881 

O.S2272 

0.B74S9 

0.09375 

0.92931 

0.97575 

0.96176 

0.95108 

0.01270 

0.02043 

0.02699 

0.03476 

0.03765 

0.03761 

0.03312 

0.03332 

0.03429 

0.03497 

0.01631 

O.OOS62 

-0.01114 

5 

39. 186 

39.602 

39.594 

39.661 

39.762 

39.836 

39.824 

39.908 

39.942 

40.003 

40.101 

40.023 

39.943 

0.09970 

0.25200 

0.35353 

0.49807 

0.63705 

0.77238 

0.81441 

0.86691 

0.89450 

0.93051 

0.97303 

0.85581 

0.93606 

0.01594 

0.02575 

0.03470 

0.04549 

0.04367 

0.04091 

0.0529S 

0.04472 

0.04520 

0.04422 

0.00300 

0.00381 

0.00558 

6 

49.141 

49.570 

49.541 

49.837 

49.714 

49.809 

49,732 

49.891 

49.939 

60.024 

60.097 

60.030 

49.039 

0.09979 

0.25255 

0.35448 

0.49767 

0.63400 

0.76630 

0.80'654 

0.85943 

0.88981 

0.9330S 

0.95041 

0.93031 

0.93368 

0.01886 

0.03062 

0.04200 

0.05606 

0.06161 

0.0S21S 

0.03704 

0.05605 

O.OSS39 

0.05127 

-0.00503 

0.01202 

0.02442 

7 

59.121 

59.554 

59.502 

59.536 

69.668 

59.777 

59.740 

59.874 

59.940 

60.052 

60.070 

60.039 

59.963 

0.09965 

0.25315 

0.35539 

0.49729 

0.63104 

0.76060 

0.70910 

0.85210 

0.00557 

0.93681 

0.93792 

0.94S40 

0.94085 

0.02 1S4 

0.03311 

0.04879 

0.06605 

0.07359 

0.07442 

0.08147 

0.00735 

0.08491 

0.0S5S9 

•0.00290 

0.01480 

0.03SS8 

a 

69. 115 

69.551 

69.481 

69.468 

69.623 

69.754 

69.702 

69.861 

69.952 

70.032 

70.059 

70.046 

09.073 

0.09994 

0.25330 

0.35636 

0.40682 

0.62819 

0.7S526 

0.79222 

0.84529 

0.80202 

0.93856 

0.92211 

0.94027 

0.94701 

0.02403 

0.03921 

0.05477 

0.07787 

0.08572 

0.08650 

0.09614 

0.07904 

0.07392 

0.05847 

0.00320 

0.01730 

0.03807 


PSI 

R 

Z 


o 9 

-n 3t 

-o? 

Q > 
?a r 

O 
^ r 

S' 


I 

G9.602 

0.43720 

0.06098 


69.760 70.050 
0.81049 0.93590 
0.09036 0.01840 


S9.520 

0.43303 

0.07532 

89.694 
0.79760 
0. 1 1895 

30.056 

0.92485 

0.02570 


3 

109.444 

0.44032 

0.03938 

109.660 
0.78718 
0. 1484S 

110.072 

0,91416 

0.03339 


4 

129.350 
0.44268 
0. 10304 

129.658 
0.77375 
0. 17848 

130. 1 14 
0.90374 
0.04153 


S 

149.209 
0.44459 
0. t 1636 

149.663 

0.77222 

0.20860 

ISO. 187 
0.89352 
0.04987 


6 

168.963 
0.44639 
0. 12953 

169.621 

0.76754 

0.23887 

170.281 

0.88370 

0.05741 


7 

188.393 
0.44710 
0. 14491 

169.413 

0.76425 

0.27033 

190.262 

0.87492 

0.07239 

PSl-180: Z(TIP)«0.061 

R(TIP)«0.879 

6 

207.821 
0.44672 
0. 16268 

209. 187 
0,76 190 
0. 30292 

210. 184 
0.66721 
0.09449 

■ 


9 


227.568 

229.043 

330, 179 

0.44618 

0,76019 

0.66001 

0. 1G04S 

0.33470 

0.11430 

10 

247.418 

248.940 

250. 166 

0.44528 

0.75872 

0.85346 

0. 19705 

0.36537 

0. 13454 

1 1 

267.288 

268.070 

270. 157 

0.44368 

0,75773 

0.84735 

0.21500 

0.39409 

0. 15368 

12 

287. 189 

208.825 

290. 155 

0.44 184 

0.75765 

0.8430S 

0.23179 

0.42339 

0. 17245 

13 

307. 104 

308.756 

310. 158 

0.44023 

0.75803 

0,83872 

0.248 13 

0.45126 

0. 18003 


t.% 




14 

326.999 

0.43888 

0.26419 

15 

346.825 

0.43761 

0.26040 

16 

366.533 

0.43603 

0.29741 

17 

366.215 

0.43400 

0.31532 

ta 

405.983 

0.43162 

0.33352 

19 

425.811 

0.42966 

0.35160 

20 

445.667 

0.427SS 

0.36936 

21 

465.548 

0.42548 

0.36666 

22 

465.450 

0.42354 

0.40365 

23 

505.350 

0.42174 

0.42048 

24 

525.219 

0.41996 

0.43749 

25 

545.040 
0.4 1803 
0.45494 


328.674 

0.75821 

0.47907 


348.578 

0.75814 

0.50738 


368.455 

0.75795 

0.53676 


388.318 

0.75788 

0.56732 


408. 187 
0.75829 
0.59868 


428.069 

0.75944 

0.63029 


447.978 

0.76100 

0.65926 


467.898 

0.76164 

0.68215 


487.845 

0.76035 

0.70163 


507.003 

0.75906 

0.72158 


527.780 

0.75623 

0.74248 


547.780 

0.75233 

0.76469 


330. 166 
0.83469 
0.20873 


350. 160 
0.83099 
0.22610 


370.044 

0.82756 

0.24464 


389.926 

0.82438 

0.26463 


409.886 

0.82141 

0.28437 


429.863 

0.81866 

0.30354 


449.857 

0.81659 

0.32219 


469.864 

0.81547 

0.34040 


489.853 

0.81434 

0.35824 


509.830 

0.81418 

0.37579 


529.791 

0.81342 

0.39341 


549,713 

0.81253 

0.41172 


PSI-360: Z(TIP)-0.235 

R(TIP}«0.829 


PS1-B40: Z(TIP>0.402 

R(TIP)-0.813 
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26 

564.640 567.731 569.637 

0.41592 0.74754 0.81142 

0.47290 0.7C345 0.43077 

27 

584.655 587.789 589.598 

0.41378 0.74225 0.81007 

0.49115 0.81364 0.44999 

28 

604.490 607.760 609.574 

0.41166 0.73691 0.60851 

0.50341 0.83873 0.46907 

29 

624.333 627.704 629.571 

0.40924 0.73095 0.80808 

0.52732 0.86195 0.48780 

30 

644.182 647.640 649.597 

0.40613 0.72367 0.81007 

0.54453 0.88270 0.50595 

31 

664.087 667.653 669.560 

0.40263 0.71605 0.81327 

0.56120 0.90171 0.52358 

32 

684.031 687.708 689.498 

0.39916 0.70918 0.81628 

0.57767 0.91969 0.54096 

33 

703.982 707.775 709.421 

0.39579 0.70310 0.81902 

0.59418 0.93690 0.55837 

PSI*‘?20: 2(T1P)»O.507 

34 R(TIP)>0.820 

723.915 727.828 729.329 

0.39247 0.69773 0.82147 

0.61089 0.95347 0.57609 

35 

743.825 747.850 749.248 

0.38922 0.69300 0.82365 

0.62784 0.96945 0.59412 

36 

763.719 767.848 769.194 

0.33614 0.68882 0.82553 

0.64494 0.98481 0.61215 

37 

783.607 787.826 789.160 

0.38334 0.68503 0.82708 

0.66204 0.99948 0.62993 
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